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ABSTRACT 

We report the discovery of two iron-group enhanced high-metallicity Mgii absorbers in 
a search through 28 Lya forest clouds along the PKS 0454+039 sight line. Based upon our 
survey and the measured redshift number densities of Wr(Mgii) < 0.3 A absorbers and Lya 
absorbers at 2; ~ 1, we suggest that roughly 5% of Lya absorbers at z < 1 will exhibit "weak" 
Mgll absorption to a ha H^i.(A2796) detection hmit of 0.02 A. The two discovered absorbers, 
at redshifts z — 0.6248 and z — 0.9315, have Wi-(Lya) = 0.33 and 0.15 A, respectively. 
Based upon photoionization modeling, the Hi column densities are inferred to be in the range 
15.8 < logA^(Hi) < 16.8 cm^^. For the z — 0.6428 absorber, if the abundance pattern is solar, 
then the cloud has [Fe/H] > —1; if its gas-phase abundance follows that of depleted clouds in our 
Galaxy, then [Fe/H] > is inferred. For the z = 0.9315 absorber, the metallicity is [Fe/H] > 0, 
whether the abundance pattern is solar or suffers depletion. Imaging and spectroscopic studies 
of the PKS 0454+039 field reveal no candidate luminous objects at these redshifts. We discuss 
the possibility that these Mgii absorbers may arise in the class of "giant" low surface brightness 
galaxies, which have [Fe/H] > —1, and even [Fe/H] > 0, in their extended disks. We tentatively 
suggest that a substantial fraction of these "weak" Mgil absorbers may select low surface 
brightness galaxies out to z ^ 1. 

Subject headings: galaxies: interstellar medium — galaxies: evolution — quasars: absorption 
lines 
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Introduction 



MetaHine absorption in the intergalactic medium, or IGM,^ is astrophysically interesting because the 
absorption properties can be exploited to reveal the star formation, chemical enrichment, and ionization 
histories of the universe. This provides a motivation for studying metal lines in Lya absorbers over as wide 
a redshift range as possible and for sampling transitions covering as many ionization levels and chemical 



species as possible (cf. Hcllstcn et al. 1997 ; Ranch, Hcahnelt, fc Stcinmetz 1997 ). Observations ( Hu et al 



1995 



1997 



Lu et al. 1996|; Kim et al. 1997) and numerical simulations ( Miralda-Escude et al. 1996; Zhang et al 



Dave et al. 1997; Norman et al. 1997) have revealed that the forest is rapidly evolving with redshift 



from z ~ 4 to z ~ 1, and that the absorbing gas is housed in a wide range of cosmic structures undergoing 
a wide range of dynamical processes. At z ~ 2, Lya clouds contain the majority of the baryon content of 
the universe. At lower redshifts, Lya clouds are thought to be more directly associated with low surface 



brightness and/or dwarf galaxies ( Salpeter 1993 ; 3hull et. al. 1996 ; Linder 1997 ), with the outer disks and 
halos of high surface brightness galaxies ( Lanzetta et al. 1995 ; Lc Brun et al. 1996|) , or with the remnant 

199^ ; 



material left over from the formation of galaxies and/or small galaxy groups (van Gorkom et. al. 
Bowen, Blades, & Pettini 1996; Lc Brun et al. 1996). Studies of the metal content and ionization conditions 



in these low redshift forest clouds, especially in the context of their association (or lack of association) with 
galaxies, could provide the "missing-link" evidence necessary for inferring the evolving interplay between 
the IGM and galaxies or the presence of low surface brightness galaxies at higher redshifts. 

A limited number of strong metal-line species have now been seen in high ionization transitions 



at z ^ 2.5 (Tytler et al. 1995; Cowie et al. 1995; Songaila & Cowie 1996). However, the chemical and 
ionization conditions of Lya clouds at low redshifts (z < 1) remain unexplored because they require 
time-intensive programs using HST. Relative to z ~ 2.5, the meta-galactic UV background flux (UVB) at 
z < 1 is reduced by a factor of '--^ 5 and its shape may be softened by stellar photons escaping bright field 



galaxies (Dcharveng et al. 1997; Giallongo, Fontana, & Madau 1997; Bergeron et al. 1994, and references 



therein). Thus, the IGM ionization conditions may have evolved so that low ionization species, especially 
the resonant Mgii AA2976, 2803 doublet and several of the stronger Fell transitions, are detectable in Lya 
clouds. As discussed below, these particular species are well suited for understanding chemical enrichment 
histories. 

Songaila & Cowie (1996, hereafter SC96) detected Civ absorption in ~ 75% of all Lya clouds at 
z ~ 2.5 with logiV(Hi) > 14.5 cm~^ and concluded that roughly 50% of log A''(Hl) < 14.5 cm~^ clouds 
could have primordial abundances. They also reported Siiv and Cii absorption in a fraction of the Lya 
clouds (including "partial" Lyman limit systems). Based upon the photoionization models of Bergeron 
& Stasihska (1986), SC96 find the metallicity of z ~ 2.5 Lya clouds to be [Z/Zq] ~ —2 and to be fairly 
uniform, with about 1 dex of scatter.^ They also report [Si/C] ratios consistent with Galactic Halo stars 
(metal poor late-type stars), in that the a-group silicon is enhanced by a factor of three over the carbon. 
This conclusion, however, is sensitive to the assumed UVB continuum shape, especially the question of how 
much star bursting galaxies contribute to the UVB, and its non-uniformity, at higher redshifts ( |Giroux"& 
ShuU 1997| ). 



■^Throughout this paper, we use the terms "Lyo cloud", "forest cloud", and "IGM" somewhat interchangeably to designate 
Lya absorption with T912 < 1. 



^Throughout this paper, we use the notation [Z/ZQ] = log Z - 
Y are any two elements. 



;Z0, and [X/Y] = log(X/Y) - log(X/y)0, where X and 
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Considering the mechanisms and range of environments that could plausibly give rise to metals in 
what are traditionally known as Lya clouds, it is difficult to understand a high level of uniformity in 
their chemical enrichment histories. As proposed by Tytler et al. (1995), there are at least three obvious 
mechanisms for the enrichment. 

(1) The larger A'^(Hi) clouds may be gravitationally bound with internal gravitational instabilities in 
which they produce their own stars, which in turn distribute the metals throughout the cloud. This type of 
object has little distinction from a galaxy. This in situ process would likely give rise to a strong metallicity 
dependency with Hi column density, unless a well-tuned mechanism governing star formation yielded a 
uniform chemical enrichment history of the IGM, as suggested by Cowie et al. (1995). Such a mechanism 
would likely represent non-standard star formation processes. 

(2) The metals may be produced in protogalaxies and then be widely distributed via mechanical ejection 



from merging events (Gnedin & Ostriker 1997) or from correlated supernovae (SNe) (Cen & Ostriker 1992). 
This implies that Lya clouds formed after the metals were distributed around the metal producing galaxies. 
The scenario also predicts that the metal enriched Lya clouds, as opposed to "Lya-only" clouds, would 
cluster like galaxies. 

(3) Population III stars, formed at z > 10 and somewhat uniformly spread throughout the IGM, may 
have distributed metals into the IGM prior to the first protogalaxies. A population of Lya forest clouds 
that have been enriched by Population III stars may exhibit IGM chemical conditions that are relatively 
unchanged from the epoch of the first stars. If so, this population would be ideal for studying the intensity 
and continuum shape evolution of the UVB from < z < 4, since the changing ionization conditions could 
be used to deconvolve the non-evolving chemical conditions from the evolving UVB. As such, the detection 
of extremely metal poor stars in the Galaxy halo would also be very interesting, since their presence would 
suggest the presence of Population III stars (cf. Ostriker fc Gnedin 1996| ). 



The chemical abundance pattern can serve as a clue to the origin, physical environment, and chemical 
enrichment history of any given Lya absorber. There are at least two major uncertainties involved in 
measuring relative abundances using QSO absorption lines: ionization corrections and dust grain depletion 
patterns. The ionization correction provides good reason for observing a wide range of ionization levels. 
Moreover, the ionization corrections are sensitive to the intensity and continuum shape of the ionizing 
radiation, possibly providing even further leverage for understanding local environments and chemical 
enrichment history. Dust depletion does not effect all a-group elements (for example, sulfur is not readily 
incorporated onto dust grains) nor all Fe-group elements (for example, zinc). However, those elements 
provide neither the strong UV absorption lines needed to accurately probe clouds with T912 < 1, nor the 
transitions observable from the ground for redshifts below z ^ 1. The strongest observable transitions are 
the Mgii AA2976,2803 doublet (a-group), and Fell A2344, 2382, and 2600. Unfortunately, both magnesium 
and iron can deplete onto dust grains and their depletion levels are environment dependent. 

For each of the three scenarios suggested above, it is expected that the relative elemental abundances 
of an enriched cloud should reflect the a-group enhanced yield of Type II SNe (note that this is consistent 
with the results of SC96). In essence, the picture is simple: if the a-group elements are enhanced relative to 
the Fe-group then the chemical enrichment has been dominated by Type II SNe. Based upon the [Si/Fe], 



[S/Fe], [0/Fe] ratios measured in Galactic halo stars, this pattern is seen for [Fe/H] < — 1 (see Lauroesch 



|et al. 1996 ). In the case of scenario (3) presented above, it is likely that only a single burst, or episode, of 



star formation would have occurred and that the metal production arise exclusively from Type II SNe. In 
the case of scenario (2), the Lya clouds would be far from the galaxies; the only metal enriched gas that 
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could infiltrate clouds forming from the primordial IGM would necessarily be ejected from correlated Type 
II SNe bursts. 

If the abundance pattern is more in line with solar proportions (i.e, [Mg/Fe] ~ [Si/Fe] ~ [Fe/H] ~ 0), 
then the picture is that Fe-group elements have been built up over a longer time scale via Type la SNe. 
This implies a star formation history local to the cloud that would have been relatively quiescent for ~ Gyr 
prior to the epoch of the observed absorption. Thus, if a given Lya cloud is measured to have [Fe/H] > — 1 
and [a/Fe]~group abundance ratios near solar proportions, then one might infer that Type la SNe have 
played a role in the cloud's chemical enrichment history. However, based upon uncertainties in Type II SNe 
yields, Gibson, Loewenstein, & Mushotzky (1997) have cautioned that the relative importance of Type la 
and Type II SNe as inter-cluster polluters remains uncertain. 

A type of extended gas-rich object that is seen to have [Fc/H] > —1, and even [Fc/H] > 0, is the 



class of giant low surface brightness galaxies (Bothun, Impey, & McGaugh 1997; Pickering & Impey 1995 



McGaugh 1994). At low redshifts, the general population of low surface brightness galaxies are seen to 



outnumber high surface brightness galaxies by a factor of at least two (Dalcanton ct al. 1997). As such, 
these galaxies are important tracers of low density dark matter halos and structure formation from small 
over-density fluctuations. They also may represent environments where the pathways of star formation and 
chemical evolution reflect non-standard astrophysical processes (Bothun et al. 1997). If these objects are 
observable in absorption over a wide range of redshifts, they likely will provide us a unique astrophysical 
laboratory for broadening our present perspective on star and galaxy formation. 

In this paper, we report the search for and discovery of Fell and Mgii doublet absorption in the 
Lya forest along PKS 0454+039 sight line. In §||, we describe the data and analysis. In we describe 
the sample of Lya lines and our search method. The absorption properties of the detected systems are 
presented in §H We apply photoionization models to the detected metal-line systems in §||, and briefly 
discuss the model results in §^ The implications of the results are addressed in §0. A brief summary is 
provided in §@. 



2. Observations and Data Analysis 



The Mgii and Fell transitions were searched for in an R = 45, 000 optical HIRES ( |Vogt et al. 1994D 



spectrum. The Lya line list was obtained from the R = 1300 G190H and G270H FOS/HST spectra of 
Boisse et al. (1997, hereafter BBLD). Three images of the PKS 0454+039 field have been incorporated into 
our study so that we may attempt to identify the luminous objects giving rise to the absorption. Two are 
high-spatial resolution WFPC2/HST images using the F450W and F702W filters (Le Brun et al. 1997, 
hereafter LBBD). The third is the deep i?[6930/1500] image (centered on A6930 with a FWHM band pass 
of 1500 A) taken from Steidel et al. (1995). We also draw upon the published (Bteidel et al. 1995) and 



unpublished (C. Stcidcl, private communication) spectroscopic redshift measurements of the many objects 



along the line of sight to the QSO. 

The HIRES spectrum was obtained and reduced as described in Churchill (1997a) and in Churchill, 
Vogt, & Charlton (1998). The HIRES spectrum has wavelength coverage 3767— 6198 A, though there are 
some breaks in the coverage redward of 5100 A because the single setting of the 2048X2048 CCD did not 
capture the complete free spectral range at these wavelengths. The FOS spectra were obtained, reduced, 
and the list of Lya forest lines used for this study were produced as described in BBLD. The acquisition and 
analysis of the high-spatial resolution WFPC2 images are described in LBBD. They also present a synopsis 
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of candidate galaxies along the line of sight to the quasar available from their study and the literature. The 
i?-band ground-based image and the spectroscopic identifications of objects in the PKS 0454+039 field are 
described in Steidel et al. (1995). 



3. Searching the Forest 

We searched the HIRES/Keck spectrum of PKS 0454+039 for Mgii absorption in the Lya 
lines reported by BBLD. In the FOS spectra, the redshift range over which Lya was detected is 
0.4163 < z(A1216) < 1.3431. The redshift range over which the Lya A1215 transition could have been 
detected was 0.41 < z < 1.69. The line list is presented in Table 0. In all cases, they are Lya-only 
systems (no other corroborating transitions are detected in absorption) . We have included only those lines 
from BBLD that arc not confused by blending or are not coincident in wavelength with strong metaHine 
transitions from the four known metal-line systems {z = 0.072,0.859, 1.068, 1.153) along the line of sight. 
The first two columns of Table Q are the redshift of Lya absorption and the rest-frame A1215.67 equivalent 
width, respectively. In columns 3-5, values of the neutral hydrogen column density, log A^(Hi), are tabulated 
for Doppler b values of 80, 30, and 15 km s~^, respectively. These are shown only to illustrate the plausible 
A^(Hi) range that might be inferred from the equivalent widths. A 6 = 30 km s^^ is representative of the 
median b value of 34 km s^^ found by Kim et al. (1997) at redshifts 2 < z < 3. There is evidence that the 
median 6 value increases with decreasing redshift, which is believed to be due to kinematic substructure 
evolution for 7V(Hi) < 14.0 cm~^ clouds. For the higher column density clouds, the widths likely reflect 



non-shock heated cloud temperatures (Hachnelt, Steinmctz, fc Rauch 1996). A 6 = 80 km s ^ is plausible 



for kinematically broadened and/or shock heated clouds ( Kim et al. 1997 ). Doppler parameters greater 



than ~ 80 km s ^ are likely due to blending (Lu ct al. 1996), possibly of physically distinct systems. The 



15 km s ^ value is the "lower cut off" value found at high redshift (Lu et al. 1996, but also see Hi; 



et al. 1995) 



The sensitivity of the search, as a function of redshift, is quantified in terms of the rest-frame limiting 
equivalent width of the Mgii A2796 transition. In Figure |l], we present the sensitivity curve, where we 
have chosen to use a 5a significance level. The redshift range over which Mgii doublets could be detected 
in the HIRES spectrum is 0.3466 < z(A2796) < 1.2134. There are small gaps in the coverage above 
2;(A2796) = 0.8340 that increase with increasing redshift. The 5a observed equivalent width detection limit 
ranged from 0.007 to 0.020 A, except for z(A2796) < 0.4662, where it ranges from 0.020 to 0.035 A. The 
results, including the detection limits and the limiting column densities of A4gii, which are based upon 
linear curve of growth analysis, are presented in columns 6 and 7 of Table 0. Only those redshifts for which 
both transitions of the Mgii doublet could be observed are tabulated. 

Absorption features were defined using our own interactive software, which is based upon the detection 
algorithms of the QSO Absorption Line Key Project ( [Schneider et al. 1993 ). The criteria that define a 



confirmed Mgii doublet are presented by Churchill et al. (1997), who have searched 26 HIRES/Keck QSO 
spectra for weak Mgii systems. These lines are fit with Gaussians to obtain their equivalent widths and 
observed central wavelengths. To locate candidate Mgii doublets, the candidate A2796 line centroid and 
detection aperture (full width at the continuum) is shifted to the expected location of the A2803 line and the 
significance level is measured. An acceptable candidate for the weaker member of the doublet occurs when 
the detection significance level is greater than or equal to the ratio of the transition /A times the significance 
level of the stronger member. Then, a "chance probability" is computed by scanning the spectrum with the 
detection aperture over the spectnmr for 50 A to both sides of the candidate and computing the fraction 
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of detected features (both emission and absorption) with a significance level greater than or equal to the 
candidate A2803 line. Most bonafide Mgii doublets have chance probabilities of < 10~^, though a very few 
have probabilities as large as ~ 10"'^. 

To a 5(7 limit of log-/V(Mgii) ^ 11.3 cm~^, the Mgii doublet was detected in two of the twenty-eight 
Lya lines in the list, which is a success rate of ~ 7%. The two weak Mgii systems found in the HIRES 
spectra, have z = 0.6428 and z = 0.9315. The data for these systems are presented in Figure ^ and their 
measured properties are listed in Tables I and |. For both, the Mgii doublet and at least one transition of 
Fell was detected. Below, we describe the measured properties of the two detected systems. 



4. Properties of the Absorbers 

In this work, we concentrate on the properties of the two absorbers for which Mgii has been detected. 
Here we note that the Lya equivalent widths are among the smallest in the sample of 28 (there are seven 
as small as or smaller than the z = 0.6428 absorber and two as small as or smaller than the z = 0.9315 
absorber). Also, we note that the Mgii and Fell transitions have been detected a factor of five to ten above 
the 5cr detection limits of the HIRES spectra. Given the stringent limits on the Mgii column densities 
for the remaining Lya absorbers, and the fact that majority appear to have higher A^(Hi) than the two 
exhibiting Mgii, it may be that there is a large dynamic range in the A(Mgii)/A^(Hi) ratio. As noted in 
, this is not the case for Civ absorption in the Lya forest. 

However, we note that upper limits on the mctallicities are not very restrictive if we assume a typical 
h parameter of 30 km s^^ and photoionization by the UVB (using CLOUDY; see Figure 11 of Churchill 
et al. 1997). For a Lya cloud with Wr(Lya) ~ 0.6 A and A^(Mgii) < 11.4 cm-^, the upper limit on [Z/Zq] 
ranges from —0.5 to —2.8, depending upon ionization level. For a Lya cloud with Wi(Lya) 0.4 and 
iV(Mgii) < 11.4 cm~^, the upper limit ranges from +0.4 to —2.2 as the cloud becomes highly ionized. 
Thus, not a great deal can be said about the range of mctallicities in the Lya forest at z < 1 based upon 
our Mgii upper limits. 



4.1. The z = 0.6428 System 

In the left hand panels of Figure ^, the Mgii and Fell HIRES profiles are presented. Also shown (top 
panel) is the FOS spectrum, with the position of the corresponding Lya line marked with a tick. Along with 
the Mgii doublet. Fell A2383 and A2600 were clearly detected. The weaker Fell transitions were covered by 
the spectrum, but were not found to the 5ct significance level. In part, this is due to the decreasing signal to 
noise below 4000 A where the HIRES sensitivity drops rapidly (see Figure |l|) . In the upper panel of Table ^, 
the measured rest-frame equivalent widths, column densities and Doppler h parameters are presented. 

The Mgii doublet ratio is 1.4 ± 0.1, and both the Mgii and Fell lines are partially resolved. Based 



upon the apparent optical depth profiles (cf. [Savage fc Sembach 199l| ), there is evidence for unresolved 



saturation in the Mgii doublet. It may be that there are two or more very narrow absorbing components 
giving rise the profile, but the signal-to-noise ratio is not high enough to model the data to this level. 
The column densities and h parameters are obtained using Voigt Profile (VP) fits that incorporated both 
the atomic physics and the instrumental spread function. We used the program minfit ( Churchill 1997a[ ), 



which performs an iterative minimization between the data and the model spectra [see Churchill (1997a) 
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for a detailed description of the convergence criteria and the error computations]. A VP model with two 
components was handed to minfit, but it returned a single component model based upon the criterion that 
there was no statistically significant difference between the best fit single component and two component 
models. The measured VP column densities for the HIRES profiles are logiV(Mgii) — 12.74 ± 0.02 cm~^ 
and log A^(Feli) = 12.46 ± 0.06 cm~^. Their respective b parameters are 6tot(Mgii) = 5.7 ± 0.3 km s~^ and 
&tot(Fell) = 4.3 ± 1.3 km s~^. The reduced xt for the simultaneous VP fit to the Mgii doublet and the two 
Fell transitions is 0.96, where the degrees of freedom is v = 107. In principle, the contribution of turbulent 
broadening to the profiles, 6turb/^tot, could be determined from the ratio of the atomic masses of iron and 
magnesium (see eq. [|l| and eq. [^), but the uncertainties are too large to directly place useful limits on 

^turb- 

As we will discuss below, well determined uncertainties in the VP quantities are central to constraining 
the ionization, thermal, and chemical conditions in the absorbing gas "cloud" . VP fits are particularly 
robust for profiles in this regime of column density and width. As shown in Churchill (1997a), the quoted 
uncertainties in the VP quantities are consistent with the spread in these quantities measured from VP fits 
to 1000 simulated spectra with similar signal to noise ratios. Thus, the measured column density and line 
broadening, and their uncertainties, are considered to be robust. 

In the FOS spectrum, the Lya absorption is the central line in a triple-blend feature. In Figure ^ 
we show the deblending fit. The adjacent lines in the blended feature are Ly7 at z = 1.1536, and Lya 
at z — 0.6448. The lone feature at 2004.5 A is Ly6 at z = 1.1536, which is the redshift of a strong Mgii 
absorber studied elsewhere (Churchill et al. 1998). Though the fit is not unique, we have adopted the 
presented result, which yielded a rest-frame Lya equivalent width of Wr = 0.33 ± 0.03 A. Because of 
the low resolution of the FOS spectrum, we could not obtain estimates of the Hi column density and b 
parameter directly from the data. There is no coverage at the expected position of the Lyman limit, so 
we cannot place an observed upper limit on A'^(Hi). The Ly/3 falls just blueward of the Lyman limit at 
1695 A due to the z — 0.8596 damped Lya absorber, and thus cannot be detected. If the Hi and the 
Mgll arise in the same physical locations in the absorber, then the total Hi b parameter is constrained to 
be 5.7 < 6(Hi) < 29.8 km s^^ (see Figure ^) including the spread introduced by the uncertainty in the 
measured 6tot(Mgii). The lower limit corresponds to the case in which turbulence or bulk motions dominate 
the line broadening and the upper limit corresponds to a thermal motions scaling, -\/246tot(Mgii). From the 
curve of growth, we have estimated that the inferred 6(Hi) range translates to a neutral hydrogen column 
density range of 14.2 < logA^(Hi) < 17.6 cm~^. For this estimate, we have included the spread introduced 
by the uncertainty in the measured IVr(Lya), which dominates over the uncertainty in 6tot(Hi). 

The FOS spectrum also covers several other transitions from a variety of species and over a wide range 
of ionization potentials. It is important to thoroughly check for the presence of absorption from these 
species and to place limits on their column densities when no absorption is detected. These limits may be 
useful for further constraining the chemical and ionization conditions of the absorber, even if the sensitivity 
level of the FOS spectrum is not very high. Thus, we have systematically searched the FOS spectrum 
for other transitions associated with the z — 0.6428 absorber, using the detection method described by 
Schneider et al. (1993). Details of the search are presented in Appendix ^ and selected results are tabulated 
in Table |. 
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4.2. The z = 0.9315 System 

In the right hand panels of Figure the Mgii and Fell HIRES profiles are presented. Also shown (top 
panel) is the FOS spectrum, with the position of the corresponding Lya line marked with a tick. Along 
with the Mgii doublet, the Fell A2383 transition was detected. The Fell A2600 transition may have been 
measurable as well, but its predicted location coincided by chance with that of the pen mark (i.e. "The 
Blob") on the HIRES Tektronic's CCD. The weaker Fell transitions were covered by the spectrum, but 
were not found to the 5(t significance level. Their limits are consistent with the Fell A2383 detection. 

In the lower panel of Table the measured rest-frame equivalent widths, column densities and Doppler 
h parameters are presented. The measured Mgii doublet ratio is 2.0 ± 0.1, though it may be closer to 1.8, 
as we discuss below. The Mgii and the Fell lines are unresolved. For unresolved lines, it is difficult to 
accurately determine the column densities and h parameters from profile fitting. 

We have performed extensive VP fitting simulations of the Mgii doublet for this system. The 
constraints for adopting the best model were the measured Mgii A2796 equivalent width and the doublet 
ratio. Using the curve of growth, we explored a grid of column densities and h parameters that were 
consistent with the measured equivalent width of the Mgii A2796 transition. The grid range was 
12.1 < log A^(Mgii) < 14.2 cm~^, corresponding to 4.0 > 6(Mgii) > 0.4 km s^^. The increments in column 
density were 0.1 dex. For each grid location, 500 spectra were simulated, convolved with the HIRES 
instrument spread function, sampled at the HIRES pixelization, and degraded to the signal-to-noise ratio 
of the observed data. The simulation output consisted of the Mgii A2796 equivalent width, the doublet 
ratio, and the VP column densities and h parameters from minfit. 

For 5 < 1 km s~^, we found that we could not recover the measured equivalent width, nor the doublet 
ratio; they both decrease dramatically with decreasing h. This is due to the finite pixelization of HIRES. As 
h is reduced, the line depth increases. As the pre-instrumentally broadened line width drops below that of 
a single pixel, saturation losses dominate. It could be argued that the measured equivalent width already 
reflects this (that the absorption is actually stronger than the measured value) and that larger equivalent 
widths should be explored. Fortunately, the VP fits recovered the input 6 values to both high accuracy 
and precision over the full range explored. For 6 < 1 km s~^, no matter the value of the equivalent width, 
the doublet ratio could never be made consistent with the data (within 3cr). For the lower limit on 6, we 
adopted the criterion that the measured doublet ratio would be a 3cr outlier in the distribution of simulated 
VP fits. For the upper limit on 6, we adopted the criterion that the measured h would be a 3(t outlier from 
the mode of the fitted h distribution. 

A caveat is worth noting. We also explored simulations in which the Mgii doublets were fit individually, 
rather than simultaneously From this wc concluded that the observed A2803 transition is likely compromised 
by a possible fiat fielding artifact in its blue wing. This is consistent with the visual appearance of the 
data in comparison to the many simulated spectra and with the measured for the VP fits. The value 
of was 1.29 with = 74 when all three transitions were fit simultaneously. If just the Mgii doublet 
was fit, then ~ 1.47 with v = 51. This value is dominated by the "poorer" fit to the A2803 transition, 
which by itself was Xy = 2.13 with v = 24. In contrast, the VP fit to just the A2796 transition yielded 
xl = 0.97 with V — 24. The effect of this residual flux in the A2803 transition was to push the measured b 
value down to 1.5 km s~^, when the doublet was fit simultaneously. When the observed A2803 transition 
was omitted from the VP fit to the data, the resulting Mgii b parameter was consistent with the mean of 
the simulations. Based upon these considerations, we have omitted the observed A2803 transition from the 
VP fit results; the measured b parameter used for interpreting the simulations was obtained by a fit to the 



- 9 - 



A2796 transition only. Since we have adopted the assumption that the observed A2803 transition has been 
compromised, we have adopted the "best" doublet ratio from the simulations, dr = 1.8 ± 0.1. This implies 
that the equivalent width of the A2803 line is slightly larger than that formally measured from the data. If 
the Mgii DR is 1.8, then rest frame A2803 equivalent width is ^ 0.023 A. 

The adopted VP column densities are logiV(Mgii) = 12.24 ± 0.09 cm"^ ^nd logiV(Feii) = 
12.29 ± 0.08 cm~^. Their respective b parameters &tot(Mgii) — 2.2 ± 0.5 km s~^ and 
&tot(Feii) = 2.3 ± 1.6 km s~^. As with the z — 0.6428, the value of 6tmb/6tot for the system could 
be determined from the ratio of the atomic masses of iron and magnesium (see eq. Q and eq. |^ ) , but the 
uncertainties are too large to directly place useful limits on btuib- 

Because of the low resolution of the FOS spectrum, we could not obtain estimates of the Hi column 
density and b parameter directly from the data. However, the wavelength range over which the Lyman 
limit break could be observed is present in the spectrum at 1760.9 A. There is no apparent flux decrement 
at the expected position of the break. However, the signal-to-noise ratio is low, ~ 5, and this places a 3a 
limit of 1.6 on the flux ratio across the break. This corresponds to an upper limit logA^(Hi) ~ 16.5 cm^^. 
The Ly/3 transition is covered at 1981.2 A, but the region is dominated by the z = 1.1537 Ly9 line at 
1983.1 A, so Ly/3 does not provide a constraint on iV(Hi). If the Hi and the Mgii arise in the same physical 
locations in the absorber, then the total Hi b parameter is constrained to be 1.8 < 6(Hi) < 13.0 km s^^ 
(see Figure ^), including the spread introduced by the uncertainty in the measured 6tot(Mgii). The 
lower limit corresponds to the case in which turbulence or bulk motions dominate the line broadening 
and the upper limit corresponds to a thermal motions scaling, •\/246tot(Mgii). From the curve of growth, 
we have estimated that the inferred b range translates to a neutral hydrogen column density range of 
13.6 < logA^(Hi) < 16.5 cm~^, where we have adopted the upper limit from the Lyman limit break 
constraint. We have included the spread introduced by the uncertainty in the measured W,:(Lya) in this 
estimate, which dominates over the uncertainty in 6tot(Mgii). 

As with the z — 0.6428 absorber, we have systematically searched the FOS spectrum for other 
transitions associated with the z = 0.9315 absorber using the detection technique described by Schneider 
et al. (1993). Details of the search are presented in Appendix ^ and selected results are tabulated in 
Table |. 



5. Modeling The Absorbers 

In order to better understand the two absorbers, we have attempted to constrain their physical 
conditions, i.e. ionization and chemical conditions, non-thermal motions, and sizes. In particular, we are 
interested in the relationship between the ionizing flux, whether it is UVB or stellar/galaxy, and the inferred 
metallicity /abundance pattern. Taken together, constraints on these two quantities may reveal a great deal 
about the origin, history, and local environment of the absorbers. 



We have modeled the absorbers as single-phase photoionized clouds using CLOUDY ( Ferland 1996 ). 
The clouds were assumed to have constant density and plane-parallel geometry. A grid of models were 
produced; for each model cloud the specified physical conditions were (1) the ionizing continuum shape and 
intensity, (2) the abundance pattern of the metals, and (3) the cloud neutral column density, iV(Hi). These 
input quantities constitute the biggest uncertainties in modeling the absorbers. We used CLOUDY in 
optimize mode, in which the residuals between the model and the measured Mgii and Fell column densities 
were minimized. The two quantities allowed to vary (optimized) were (1) the metallicity of the assumed 
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metal abundance patterns, and (2) the total hydrogen density, rtu- For those ionization species for which 
column density upper limits were available, we applied the upper limits to the models. 

5.1. The Photoionizing Sources 

The two physical conditions within the absorbers that are the most telling of its formation history are 
their abundance pattern/metallicity and their photoionization source, either local stellar radiation or the 
UV background (UVB). In fact, the inferred chemical conditions are sensitive to the intensity and shape of 
the ionizing flux continuum. There are several scenarios and we address a few of the more obvious ones 
below. 



5.1.1. The UVB Scenario 

The two absorbers, whether associated with galaxies or not, may have photoionization conditions 
dominated by the UVB. To model this possibility, we have employed the UVB spectrum of Haardt & 
Madau (1996), where the intensity has been normalized at z = 0.5 and z — 1.0 for the z — 0.6428 and the 
z = 0.9315 absorbers, respectively. The Haardt & Madau UVB spectrum accounts not only for the UV flux 
emitted by QSOs and active galactic nuclei, but also for additional UV flux due to the reprocessing of soft 
X-rays (also from the QSOs and active galactic nuclei) in intervening absorbers at all redshifts. 

The UVB spectra are shown in Figure |[ where only a select range of energies is shown. Also illustrated 
are the locations of the ionization potentials of a few key ionization species. The relevant ionization 
potentials that we are studying are all just above 1 Ryd, the ionization potential of Hi. Mgii and Fell have 
ionization potentials of 1.11 and 1.19 Ryd, respectively. The Cii ionization potential is 1.79 Ryd, and for 
Cm is 3.52 Ryd. We mention Cii and Cm because they probe the Hei edge, at 1.81 Ryd, and because we 
have limits on the Cii and Cm absorption strengths. 

5.1.2. The Stellar/ Galaxy Scenario 

These particular clouds could be embedded within galaxies that are aligned with the QSO on the plane 
of the sky (zero impact parameter), or they could be in the outskirts of the galaxies, i.e. in the extended 
halo or outer disks. In the latter scenario, the radiating stars can be treated as if they are all equidistant 
from the clouds. For any stellar/galaxy scenario, the number of stars and their spectral types, metallicities, 
and distances (quantities that determine the intensity and continuum shape of the ionizing flux) must be 
consistent with known objects in the universe. 

The stellar/galactic UV flux could arise from a late-type solar-metallicity stellar population, which 
would give rise to a rapidly falling continuum with large Hi and Hei breaks. A "soft" spectrum is required 
by the observed upper limits on the high ionization species. A central question defining the scenario is: 
what level could stellar/galaxy flux be contributing to the UVB or be completely dominating the UVB? We 
have explored this question and have outlined the astrophysical principles in Appendix We constructed 
stellar/galactic CLOUDY grids that included three galactic spectral energy distribution models over a range 
of intensities and covered cases in which the stellar/galaxy flux was progressively stronger compared to the 
UVB and the case in which the UVB was locally "extinct" . 
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For the "dominant stellar-type" scenarios, in which the cloud could be near a dominant single star, 



we used Atlas stellar models ( Kurucz 1991 ). We produced a grid of optimized CLOUDY models using 
solar metallicity stars with Teff = 6,000, 10,000, 15,000, 20,000, 30,000 K, and logg = 4.4 (the solar value). 
The spectral shape is not sensitive to the surface gravity, but is quite sensitive to the metallicity and the 
effective surface temperature. The continuum falls more rapidly toward the UV for solar metallicity stars, 
so these stars have "softer" continua than what would be expected in a low metallicity early-type galaxy. 

To account for various metallicities and/or stellar populations, we also produced optimized CLOUDY 
grids using synthetic galaxy spectra. We used a 12 Gyr single-burst Worthey (1994) model with metallicity 
[Z/Z0] ~ —0.7, and a somewhat younger Worthey model with a 8 Gyr single burst with metallicity 
[Z/Zq] = —2. We also used a later-type galaxy model from Bruzual & Chariot (1993) with an exponentially 
decreasing star formation rate (SFR). This model is a 16 Gyr stellar population with 1% of the total 
star-forming mass in stars after a Gyr. These models serve to bracket a reasonable spread in galaxy spectral 
properties, given that an extreme scenario such as a star bursting galaxy can be ruled out for two reasons. 
First, a star burst spectrum would highly ionize the gas, giving rise to Siiv and Civ absorption out to a 
galactocentric distance of ~ 100 kpc (cf. eq. [2] of piroux fc ShuU 1997 ). In fact, we found that it was 



difhcult to not produce too much Siiv and CiV even with the exponential SFR model. Second, images of 
the PKS 0454+039 field, in which point spread function removal of the QSO has been performed to high 
accuracy, reveal no unidentified luminous objects with the characteristics of a star bursting galaxy to a 
limiting K magnitude of 20.5. 



5.2. The Metallicity and Abundance Pattern 

Different chemical enrichment histories and different environments can give rise to a wide variety of 
chemical and ionization conditions. Given the possible high iron to magnesium abundance ratio in these 
absorbers, it is reasonable to assume that the clouds could arise in or near galaxies. Thus, to better 
understand the origin of the absorbers, we modeled three abundance patterns that are taken from typical 
gaseous objects in galaxies. The first is the solar abundance pattern, taken from Grevesse & Anders (1989) 
and Grevesse & Noels (1993). 



In the interstellar medium (ISM), both magnesium and iron deplete onto dust grains (cf. Lauroesch 



et al. 1996: Savage fc SembachT99^ ). Since the main constraints on the CLOUDY optimization are the 



measured Mgii and Fell column densities, we also considered the effects of heating and cooling by grains 
and the dust depleted abundance patterns of two common interstellar environments. We used the Hii 
abundance depletion pattern taken from Baldwin et al. (1991), Rubin et al. (1991), and Osterbrock, Tran, 



& Vielleux (1992). For this abundance pattern we used the "large-R" grains (Baldwin et al. 1991), which 
are characterized by a more or less grey UV extinction. Like the solar abundances, the Hii pattern, which 
has [Mg/Fe] — —0.07, provides a good template for an iron-group enriched chemical evolution history. 

We have also explored the possibility that these absorbers may actually be a-group enhanced. Thus, 
we used the abundance patterns reported by Cowie & Songaila (1986) for the cold and warm phases of the 
ISM. This pattern is characterized by [Mg/Fe] = 1.23, and overall enhanced a-group elements. The dust 
grains used in these ISM models are a mixture of the graphite and silicates. 

The presence of grains also effects the cooling and heating balance, and thus the ionization balance, of 
the clouds. As mentioned, the metallicity, or the scaling factor of the input abundance pattern for elements 
heavier than helium, was allowed to vary. As the metallicity in a cloud is increased by factors of a few. 
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the cooling rates are dramatically increased and the cloud equilibrium temperatures drop significantly. 
For some clouds, the optimal metallicity was 10 to 100 times the initial input and the cloud equilibrium 
temperature in the 100 K range. 



5.3. Turbulent and/or Bulk Motions 

It is not possible to obtain a direct determination of the Hi column densities for the two Mgii absorbers. 
For each, the extreme range of plausible Hi column densities can be estimated from the measured Mgii 
b parameter for the assumption of a thermal line broadening [lower iV(Hi) limit] or turbulent broadening 
[upper A^(Hi) limit]. The range turns out to be large, 13.8 < log A^(Hi) < 17.0 cm^^, when the uncertainties 
in &(Hi) and Wi-{Lya) are considered. We have modeled clouds with logA^(Hi) — 14.50, 15.50, 15.75, 
16.00, 16.25, 16.50, 16.75, 17.00, and 17.50 cm^^. As will be shown, there were no logA^(Hi) < 15.75 cm^^ 
cloud models consistent with the data. Since, for a known equivalent width, the curve of growth provides a 
direct relationship between the Hi column density and b parameter, both quantities should be considered 
for constraining viable cloud models. Thus, we have developed a technique to further constrain the 
acceptable model cloud properties using the relationship between the turbulent and thermal line broadening 
mechanisms. 

Using the equilibrium temperature of the CLOUDY models and the &tot of Mgii, one can narrow 
parameter space to those models that are self-consistent in 6tuib/^tot: temperature, and A'^(Hi). Here, 6turb 
is the non-thermal contribution to the line broadening. Assuming the contributions are Gaussian, &tot is 
written. 



;,2 _ i2 I /,2 _ 2kT ,2 ^1 ^ 

"tot ~ "therm + "turb ~ ^ ^ "turbi 

m 



where T is the equilibrium kinetic temperature, and m is the ion mass. In reality, non-thermal motions are 
not expected to be Gaussian (see discussion in §5.5), so eq. ^ is simply a parameterization of the relative 



contributions. Defining / = &turb/&tot, the kinetic temperature of the cloud is written, 

r~ 1450 (l-/2)6L(Mgii) K. (2) 

In Figure ^, the predicted kinetic temperature for absorbers are shown as a function of /. Schematically 



superimposed are the kinetic temperature ranges of typical gaseous structures in the Galaxy (Savage & 



Bembach 199(: ; Fitzpatrick fc Spitzer 1994 ; Osterbrock 1989| ; jpitzer 197g). Also shown is the inferred range 



of Mgii absorbers, as observed using QSO absorption lines ( Churchill 1997a ; Churchill et al. 1998). This 
diagram is discussed further in §^ 

The b parameter of each species is required in order to obtain its estimated column density from the 
curve of growth. Usage of the curve of growth method is consistent with our assumption of single-phase 
clouds, given that the velocity structure of the low ionization absorption profiles are not complicated and 
are suggestive of a single absorbing component. Assuming that all species are subject to the same thermal 
and non-thermal conditions, the total Doppler b parameter of any species, X, can be obtained from, 

6tot(X) ^ 1^ (1 - f) + /2| ^ &tot(Mgii) km s-\ (3) 
where A is the nucleon number (^Mg — 24). 
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5.4. Application of Constraints to Models 

In the absence of direct measurements of Af(Hi), the most acceptable cloud models are those that are 
self-consistent in that their Af(Hi), temperatures, and / are simultaneously consistent with those allowed 
by the data. With the above formalism in hand, the application for constraining a given cloud model is 
as follows: The range of acceptable / values parameterize the clouds, and give the inferred non-thermal 
component to the line broadening. The range is determined from the model cloud kinetic temperature, T, 
using a modified version of eq. 

1 f T \ 



■^(^^ ^ 52 (Mgii) 1^1450; ' 



where 6(-|-)(Mgii) — 6tot(Mgii) ± cTf,. The cloud temperature was taken as a simple average, given that the 
model clouds are constant density by definition. For clouds with iV(Hi) < 16.5 cm"^^ the temperature 
was constant as a function of depth. For the higher column density clouds, ionization and temperature 
structure was present, but by no more than 10%. Eq. ||] was then used to obtain the inferred 6tot(Hi), as 
illustrated in Figure ^ and Figure Using curve of growth analysis, the range of acceptable hydrogen 
column densities is obtained from the measured Lya equivalent width (and its uncertainty) and the range 
of acceptable 6tot(Hl). 



5.5. Caveats Regarding the Model Design 

The quoted upper limits on the ionization species covered by the FOS spectrum were estimated 
assuming a single-phase isothermal cloud. There are a priori reasons that we have adopted this assumption. 
First, if the absorbers are Lya clouds that are pressure confined by a warmer, more tenuous medium (more 
than one thermal phase giving rise to absorption), they cannot be heated by the UVB ( [Guilbert, Fabian" 



fc McCray 1983| ; Donahue fc ShuU 1991 ). On the other hand, under the assumption of photoionization by 



the UVB, a very low ionization parameter, U = Uj/nn < 10 , is required in order for the absorbing 



gas to give rise to iV(Feii) ~ 7V(Mgii) (Bergeron & Stasihska 1986; Churchill 1997a). Therefore, for 



UVB photoionization, clouds giving rise to the Mgii and comparable Fell absorption must have low 
ionization conditions. If strong Civ and Siiv were to arise in a UVB photoionized absorber exhibiting 
7V(Mgii) ~ iV(Feii), then the abundance pattern would need to be very different from solar, with a 1 — 2 
dex enhancement of iron or decrement of magnesium. Common chemical evolution paths do not give rise to 
this type of abundance pattern. 

Nonetheless, the possibility remains that the absorbers could be multiphased in their ionization 
structures (perhaps they are not pressure confined or are ionized by stellar/galactic flux), so that the 
measured Lya equivalent width may be due to a kinematic complex of clouds in which only one is giving 
rise to the observed Mgii and Fell absorption. As reported by SC96, Lya clouds at z = 2.5 are seen to 
exhibit kinematic structure in Civ absorption, with the highest velocity components more highly ionized. 
They suggest that this arises due to layered stratification in coalescing clouds. If this were the case for 
the low redshift absorbers reported here, then the upper limits provided by the FOS spectrum would not 
apply to this higher ionization phase; the limits would be lower than those quoted above on account of the 
expected larger b parameters. 

If a high ionization phase exists in either absorber, it has apparently not been detected in the FOS 
spectrum, suggesting that either its ionization conditions or abundance patterns are not giving rise to 
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strong absorption lines. For the ionization conditions to give rise to a multiphase absorber, the neutral 
hydrogen column density would need to approach the Lyman limit value of lO^'' '^ cm~^ to provide shielding 
to the low ionization phase. If there is no surrounding hot phase associated with these absorbers, it might 
be difficult to understand them in terms of pressure confined clouds; they could either be gravitationally 
bound or not be in dynamical equilibrium. Observationally, some ambiguity still remains for the detection 
of Civ and S IV in the z = 0.6428 absorber. This uncertainty could be checked by higher resolution and 
higher signal-to-noise ratio observations with STIS/HST. 

A main point here is that the physical region of the absorbers giving rise to the Mgii and Fell must 
have insignificant Civ and Siiv abundances. A second main point is that the estimated Hi column density 
in this region would be an upper limit to this low ionization phase, since some fraction of Hi would be 
arising in the high ionization phase as well. Unless the path length through a possible high ionization phase 
is far greater than that of the low ionization phase, the estimated ranges of the Hi column densities for the 
two absorbers cannot be significantly incorrect. 

There is also the possibility that the spatial alignment of the neutral hydrogen in these absorbers 
may not directly coincide with that of the ionized magnesium (ionization structure), since the ionization 



potential of Mgii and Hi are slightly different (P. Boisse, private communication). Again, this would imply 



that the Hi column densities in the region where Mgii and Fell are giving rise to absorption is smaller than 
the estimated values. 

The presence of dust can effect the modeling in two ways. First, for the model clouds using the Hii/ISM 
dust-depleted abundance pattern, there is a breakdown in the inferred metallicity when the temperatures 
are below T ~ 1000 K. The condensation temperature (the temperature at which 50% of an element is 
removed from the gas phase due to depletion) of both magnesium and iron is ^ 1300 K. This means that 
these cool clouds would have a significantly higher depletion than the warmer clouds, and this has not been 
accounted for in the CLOUDY modeling. The implication is that the optimized metallicity is in fact an 
underestimate of the gas-phase abundances, given that the depletion would be larger than that input into 
the modelj^ Second, if dust is present, it can seriously modify the intensity and shape of the incident UV 
spectrum, which would have implications for the inferred source of the radiation from the photoionization 
models. We discuss this point further in §|6.3|. 



Central to the model interpretation is the assumption that the line broadening is governed by Gaussian 



distribution functions for the particle velocities. N-body simulations with hydrodynamics (Norman et al 



1997 ; Dave et al. 1997 ; [Zhang et al. 1997 ) have shown that the absorbers are filamentary structures and 



that the concept of an absorbing "cloud" is not altogether valid. Often, the simulations result in gas which 
is collapsing towards (or expanding along) a filament; this results in hydrodynamic features and shocks 
in which the distribution function of velocities is not Gaussian. If the Mgii absorption profiles give any 
indication of the velocity distribution function, then the absorbers are not inconsistent with a Gaussian; in 
fact, very quiescent gas is suggested by the HIRES profiles. 

A final caveat is that we did not include turbulence physics in the CLOUDY models (only micro- 
turbulence and not bulk motions could have been modeled using CLOUDY). Thus, turbulence is not 
treated self consistently within the framework used to infer the absorber physical conditions. The inclusion 
of turbulence would result in an additional pressure source, which would effect the balance between the 



^These T < 1000 K model clouds are already characterized by quite large metallicities, even without accounting for increased 
condensation. 
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model cloud density and its depth. However, the model cloud densities were allowed to vary in an optimized 
fashion, and since the cloud depth is adjusted with each equilibrium calculation the pressure adjustment 
has little effect. The second consequence of modeling turbulence is that the line-center optical depths of 
absorption transitions decrease, while they increase in the line wings. However, since the majority of the 
model clouds in our grids were optically thin to neutral hydrogen, any change in the line profile shape has 
little effect on the model cloud equilibrium. 



6. Model Results 



In the final analysis, only the scenario in which the model clouds arc photoionized by the UVB was 
both astrophysically plausible and consistent with the data. Here, we focus on the results for the UVB 
scenario with the solar and Hii/ISM abundance patterns and then address to what level the stellar/galactic 
scenarios can be ruled out as viable ionizing sources. We defer discussion of the implications of the model 
results until §0. 

The optimized metallicities and densities are given for both abundance patterns in Tables ^ and ^ 
Cloud models with the a-group enhanced abundance pattern (Cowie & Songaila 1986) did not converge 
within the allowed uncertainties in the Mgii and Fell column densities. Thus, we conclude that neither 
the z = 0.6428 nor the z — 0.9315 absorber have a-group enhanced abundance patterns. This implies 



a gas-phase [Fe/H] > — 1 ( Laurocsch ct al. 1996 ; Savage fc Sembach 1996 ), and possibly iron-group 
enrichment by Type la SNe, though this remains somewhat controversial (Gibson et al. 1997). The iV(Hl) 
versus / parameter space is illustrated in Figures ^ and ^ for the data and for the model clouds. The 
adopted ranges for the absorber Hi column densities are defined by the overlap of the allowed ranges 
constrained by the data and by the CLOUDY models. Interestingly, the results indicate that the clouds 
have a substantial non-thermal line broadening (i.e. they are turbulent or are undergoing differential bulk 
motions). However, it is striking that the Mgii profiles reveal a velocity dispersion of only a few km s~^. 
This provides a counter example to the expected large b parameters if the gas was not dynamically settled 
(as found in hydrodynamic simulations of the Lya forest). The very quiet nature of these absorbers indicate 
that our assumption of a well defined temperature for a settled gas "cloud" is well founded. A possible, 
though unlikely, counter example would be if the absorbers were streaming filaments seen perpendicular to 
their elongation and streaming motion. 

Shown in Figure]^ are the predicted kinetic temperatures of the absorbers as a function of &turb/^tot- 
Temperature and turbulence limits may provide clues to the nature of the absorbing gas when the range 
of inferred properties are compared to gaseous objects typically found in galaxies. Planetary nebulae have 
T > 40,000 K ( psterbrock 1989| ; ^pitzer 1978] ). If the z = 0.6428 absorber has / < 0.4, then its inferred 
temperature is consistent with that of a planetary nebulae. Typical expansion velocities of these objects 
have line widths of ~ 20 — 30 km s^"'^ ( Psterbrock 1989| ). Since the Mgii and Fell Doppler parameters 
would also reflect the expansion velocities, it is highly unlikely that the z = 0.6428 absorber arises in 
a planetary nebula. For 0.4 < / < 0.9, the z — 0.6428 absorber temperature is consistent with that of 



individual clouds in complex Mgii systems (Churchill 1997a; Churchill et al. 1998). Both absorbers are 
consistent with the temperature range of Hii regions and warm Hi clouds, 7000 < r(Hii) < 14,000 K 
and 4000 < r(Hi) < 7000 K, respectively ( |Fitzpatrick fc Spitzer 19941 [Psterbrock 1989| ; ^pitzer 1978| ). 
For the z ~ 0.9315 absorber, the inferred 6turb/&tot would be / < 0.8, and for the z = 0.6428 absorber / 
would be confined to the narrow range 0.90 < / < 0.95. For these / values, the non-thermal broadening 
would be roug hly 6turb(Mgii) < 5 km s^^ and < 2 km s^^ for the z = 0.6428 and the z = 0.9315 
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absorber, respectively. The typical sound speeds in Hii and Hi clouds are ^ 10 km s"'^ ( Spitzer 197§| ). If 
these absorbers are Hii or Hi clouds similar to those found in the Galaxy, the line of sight non-thermal 
broadening is well below that expected for propagating disturbances in the clouds. If the absorbers are 
dominated by turbulent motions (/ > 0.98), then they must have T < 150 K. Diffuse ISM clouds have 
typical temperatures in the range 30 < T < 150 K (Spitzer 1978). In this regime, it would be more likely 
that the broadening was dominated by bulk flows rather than internal turbulence, given that the turbulent 
motion would propagate at the ~ 0.1 km s^^ sound speed typical of diffuse clouds (Spitzer 1978). 



6.1. The z = 0.6428 Absorber Properties 

Assuming the solar abundance pattern, the neutral hydrogen of the z — 0.6428 model cloud is 
in the range 16.3 < logA(Hi) < 16.8 cm^^ The range of fetmb/^tot is 0.85 < / < 0.93, which 
correspond to the kinetic temperatures 13,000 >T> 6500 K. The metallicity and model cloud density are 
-0.2 > [Z/Zq] > -0.7, and 0.01 < uh < 0.02 cm"^, for the range of A^(Hi). For the Hii abundance pattern, 
the inferred neutral hydrogen column density is slightly higher, in the range 16.7 < log7V(Hi) < 17.2 cm~^. 
The range of 6turb/6tot is 0.90 < / < 0.95, making this cloud kinetic temperature somewhat lower than the 
solar abundance model. The density and metallicity are nn ~ 0.008 cm~'^ and +0.4 > [Z/Zjjjj] > 0.0, 
for the range of Af(Hi). It could be that this cloud has very enhanced metallicity with an Hii abundance 
pattern, but this is a far reaching suggestion given the range allowed by the solar abundance pattern. Still, 
the cloud is inferred to have gas-phase [Fe/H] > — 1. 

These model clouds are difficult to understand in terms of objects typical of the Galactic disk or the 
Magellanic Clouds. For one, the typical A(Hi) observed in Galactic objects is log TV > 19.5 cm~^ ( |Savage 



fe Sembach 1996), two or more orders of magnitude greater than what is inferred for this absorber. Second, 



the typical density of T ~ 10, 000 K clouds (warm low density medium) is (nn) ~ 0.2 cm ^ ( Spitzer 1978 ) 



which is higher than allowed by the optimal models. Third, the inferred gas-phase abundances in the the 



warm and cool disk are [Fe/H] ~ — 1.1 and ^ —2.1, respectively ( Savage fc Sembach 1996 ). These fall well 
below those predicted by the models. However, the metallicity range of the solar abundance pattern model 
is consistent with [Fe/H] ~ —0.6 found for the Galactic Halo. 



6.2. The z — 0.9315 Absorber Properties 

The z = 0.9315 model cloud appears to have a high gas-phase metallicity, whether it be super solar 
or an enhanced Hii pattern. Assuming the solar abundance pattern, the neutral hydrogen is in the range 
15.8 < logA^(Hi) < 16.3 cm~^. Note that this is consistent with the upper limit of 16.5 cm^^ inferred from 
the lack of a Lyman limit break in the FOS spectrum. The range of 6tuib/^tot is 0.62 < / < 0.94, which 
correspond to the kinetic temperatures 4600 > T > 870 K. The metallicity and model cloud density are 
+0.7 > [Z/Zq] > +0.1, and 0.2 < < 0.4 cm"^, for the range of A(Hi). The optimal solar abundance 
model cloud is relatively dense with up to five times solar abundance. The inferred [Fe/H] is even greater 
for the Hii abundance pattern model. The neutral hydrogen column density is slightly higher, in the range 
16.0 < logA^(Hi) < 16.5 cm^^ (also consistent with the lack of a Lyman limit). The range of ^turb/^tot is 
0.75 < / < 0.97, which corresponds to the kinetic temperatures 3300 > T > 450 K. There is an inversion in 
the cooling curve at T ~ 2000 K, which implies that this absorber cannot be stable across the full range of 
allowed temperatures. It must either be a few thousand degrees or several hundred degrees. The density 
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and metallicity are — 0.1 cm~'^ and +1.6 > [Z/Zjjjj] > +0.7, for the respective iV(Hi). This is a 
metanicity enhancement of five to 40 times over the typical values seen in Galactic Hii regions (Baldwin 

Osterbrock et al. 1992). To date, no other intervening QSO absorption 



et al. 1991; Rubin et al. 1991 



system with such a high metallicity has been reported. 



6.3. Why the Stellar/Galaxy Scenarios Fail 

In what follows, we discuss the difficulties with the stellar/galaxy scenarios. The constraints were a 
trade off between the number of stars, or their number density, and the stellar population. The former is 
constrained by astrophysics, assuming a non-extreme stellar environment, and by the imaging data. The 
latter is constrained by the absorption line data, which have limited the UV ionizing flux to late-type stars 
and/or early-type galaxies. 

First, consider the case in which the stellar/galactic flux contributes to the UVB intensity. In order 
for a stellar /galactic contribution to modify the properties of a UVB model cloud, the stellar/galactic 
flux must exceed the UVB at > 1 Ryd, particularly in the regions 1 < hv < 1.2 Ryd (from the Hi edge 
up to and including the Fell ionization potential). As outlined in Appendix |b[ if the stellar population 
is dominated by AOIII and AOV stars, this requires ^ 10^^ stars confined to a region of space ^ 1 kpc 
in radius. This implies a stellar number density n» > 500 stars pc~^, which is five orders of magnitude 



greater than the density of AOV stars in the solar neighborhood ( Allen 1981 ). The required number density 
increases dramatically for later spectral types. If, on the other hand, the stars are dominated by early-type 
BOV (BOI,III) stars, then ~ 10,000 (1000) stars would be required in a volume of radius 1 kpc. For 
the main sequence stars, this corresponds to a number density about 100 times greater than that of the 



solar neighborhood ( Allen 1981 ). Only O stars can provide the UV fiux necessary to match the UVB at 
1 Ryd and have a number density of stars consistent with that of a typical galaxy environment. However, 
early-type stars give rise to high ionization absorption properties, especially Cm, CiV, and Siiv, so that 
the model cloud conditions are not consistent with the data. 

The 12 Gyr [Z/Zq] = —0.7 Worthey (1994) model is characterized by a steep continuum slope with 
IWogvfi, ^ —5 from 1 to 1.2 Ryd. This contimmm shape is a smooth continuation of the Hi edge, so that 
this galaxy model had to have vf^, > 10* times that of the UVB at 5500 A before affecting the model 
cloud properties. Based upon the arguments presented in Appendix ^ for the expected distribution of 



main sequence and giant stars in these galaxies (Worthey 1994), the stellar number densities would be 



extreme, > 10^ stars pc or ^ 10^^ stars within a kpc. Even under these extreme conditions, the 12 
Gyr Worthey models yielded [Fe/H] > — 1. 

The 8 Gyr Worthey model with metallicity [Z/Zq] = —2 is characterized by a continuum with a 
smaller Hi edge of Alogi//,y ~ —2 at 1 Ryd and a power law with ~ out to ^ 1.5 Ryd (this is steep 
but not nearly as steep as the 12 Gyr model). The fiux must be elevated to i/f^ > 10^ times that of the 
UVB at 5500 A before affecting the model cloud properties. As the flux, i/f^{0.17), is increased above 
^ 10^ ergs cm~^ s~^, the ratio iV(Hi)/A^(Hii) very quickly changes from unity to ~ 10~^; the model 
cloud becomes highly ionized and the limits on Cm and/or Civ are exceeded (depending upon the specific 
absorber). To adopt the idea that the 8 Gyr Worthey galaxy is contributing to the ionization conditions, we 
would be forced to accept a very narrow range of acceptable z//i,(0.17) arising from the galaxy. This range 
implies the absorbers would be embedded in the galaxy itself (zero impact parameter), with an extremely 
unrealistic number of stars. Such a scenario is ruled out. 
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The exponential SFR spectrum of Bruzual & Chariot (1993) has a continuum shape at 1 — 3 Ryd that 
is similar to the Haardt & Madau UVB spectrum. Thus, as the galactic spectrum is slowly increased over 
the UVB, the model cloud properties adjust such that the ionization parameter is held constant (the cloud 
density increases, but the metallicity and temperature do not change). At j^/y(0.17) ^ 0.1 ergs cm^^ s^^, 
there is a slow decrease in the ratio iV(Hi)/iV(Hii); the model cloud becomes progressively highly ionized. 
The dominant ionization species of magnesium and iron are then Mglil and Felii, respectively, and the Civ 
and Si IV column densities exceed the limits allowed by the data. The exponential SFR galaxy of Bruzual 
and Chariot is not a viable photoionizing source. The galaxy flux makes little modification to the inferred 
cloud properties (does not modify UVB only models) for a large range of intensities; once it does modify 
the UVB models, the cloud properties are inconsistent with the data. 

Now, consider the case in which the stellar/galactic flux is dominant. We examined this by excluding 
the Haardt & Madau UVB from the incident flux. This would require environments that are more 
extreme than those presented above, or that are shielded from the UVB but not from the stellar flux. Is 
it possible that dust extinction could be playing a role? For a dust absorption/scattering cross section 
of (t(1 Ryd) ~ 10~^^ cm~^ ( [Mathis, Rumpl, fc Nordsieck 1977 ), the total hydrogen column density 



[N(Hi) + iV(Hii) + 27V(H2)] required for Tdust 3 at 1 Ryd is logiV 22.5 cm~^ Given the upper limits 
on A^(Hi) for the absorbers, this implies highly ionized gas with A^(Hi)/iVtot(H) > 10~^. In all our model 
clouds, the fraction of molecular hydrogen never increased above /(II2) < 10^^. 

In principle, one could imagine an environment in which the absorbing gas is embedded in a late-type 
dwarf galaxy that itself is enshrouded in dust. This scenario provides the steep continuum above 1 Ryd 
while not requiring the spectrum intensity to be elevated above the UVB. For some models it was possible 
to achieve logA^(Hii) ~ 22.5 cm"^, but even with no competition from the UVB, these shielded models 
required high z//i,(0.17). Again, as outlined in Appendix ^ this implies unrealistic stellar densities. It 
is not simply a matter of the continuum slope, but also the intensity that dictates the cloud ionization 
conditions. Interestingly, these models produced very low metallicities, with [Z/Zq] ^ —3. All of the above 
stellar/galaxy models represent our failed attempt to locate a place in parameter space consistent with low 
metallicity clouds. 



7. Discussion 

Based upon the success rate of finding two Mgii absorbers out of 28 Lya absorbers along the 
PKS 0454+039 line of sight, we must conclude that the two absorbers are unique in some manner with 
respect to the Lya forest at large. We note, as mentioned in that not a great deal can be quantified 
about the range of metallicities in the Lya forest at z < 1 based upon our upper limits on the Mgii 
column densities. However, if our assumption of ~ 30 km s~^ is not applicable, as found for the two Mgii 
absorbers, then our quoted estimates in the range of metallicities would be quite wrong. In fact, given that 
the upper limits on the metallicity are super-solar for the Wr(Lya) < 0.3 A absorbers, and that it is not 
expected that many of the absorbers are super-solar, we are lead to suggest that either the b parameters 
of the M^i.(Lya) < 0.3 A absorbers are significantly smaller than 30 km s^^ or the objects do not have a 
Gaussian velocity distribution (simple curve of growth techniques are not applicable). It could be that 
the uniqueness of the two Mgii absorbers is that they are dynamically settled and have a small velocity 
dispersion. 

Other possibilities for the uniqueness of the Mgii absorbers, are: 
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(1) Their ionization conditions are different: If the chemical enrichment histories of these absorbers are 
typical of Lya forest clouds, then it might be inferred that their ionization conditions are governed by a 
local UV flux rather than the UVB. However, our models lead us to conclude that there is nothing "special" 
about the photoionization conditions of these clouds; they are best described as being photoionized by the 
UVB. There is no evidence to suggest that these clouds are coUisionally ionized. 

(2) Their chemical conditions are different: If the ionization source is no different than that ionizing 
the other 26 Lya systems we searched, then we are led to infer somewhat unique chemical enrichment 
histories for these two absorbers. If these absorbers are IGM/Lya clouds, in that they are not associated 
with galaxies, then the low metallicity results of SC96 at higher redshift do not support the notion that 
these two absorbers have undergone typical IGM chemical enrichment. In other words, these two absorbers 
are not consistent with the picture in which the IGM was enriched at high redshift by a single burst of 
Population III stars that left the Lya forest a-group enhanced^] with [Fe/H] < — 2. Thus, these absorbers 
could constitute a metal-rich minority of the IGM/Lya population (less than ^ 10%). It then becomes a 
question of understanding what environments and evolutionary histories give rise to high metal content in 
some Lya clouds. It is already known that at least some fraction of the Lya forest at z < 1 are associated 



with galaxies (Le Brun et al. 1996; Lanzetta et al. 1995; Bowen et al. 1996) 



(3) Both their ionization and chemical conditions are different: Inferring the chemical content of ionized 
absorbers requires ionization corrections that are uncertain. In the case of photoionization, the inferred 
conditions of the clouds are very sensitive to both the intensity and shape of the ionizing continuum. 
We have explored this interplay between the chemical content of the clouds and the properties of the 
UV ionizing flux, and have concluded with some certainty that the ionizing field for the two absorbers is 
constrained to have slope and intensity consistent with the UVB. This is tantamount to saying that it is 
only the chemical conditions that are inferred to be unique in the two absorbers. 

Overall, it is difficult to understand these absorbers in terms of the classic picture of Lya forest clouds. 
One problem is that their inferred Hi column densities are higher than "typical" forest clouds. A second is 
their high [Fe/H] and iron-group enhance abundance pattern. Since the metal abundance patterns of these 
absorbers are not a-group enhanced, it is implied that their environments have been influenced by Type la 



SNe yields (Lauroesch et al. 1996; Timmes, Lauroesch, fc Truran 1995; however, see Gibson et al. 1997). It 



could be that these absorbers are associated with galaxies. However, based upon imaging and spectroscopic 
studies, there are no candidate high surface brightness (HSB) objects in the PKS 0454+039 field. 

The extended luminous objects identified in the WFPC2 images (Figure 4 of LBBD) and ground-based 
image (Figure 3 of Steidel et al. 1995), have now had their redshifts spectroscopically measured using LRIS 



on the Keck I telescope. None of them have redshifts that match the two Mgii absorbers ( C. Steidel 



private communication). Object #5 in LBBD has been confirmed by Steidel and collaborators to be the 
strong Mgii absorbing galaxy at z = 1.1536. At impacts greater than 120h^^ there are three galaxies with 
0.605 < z < 0.610 (not presented in published images). Thus, to a limiting of magnitude oi K ^ 20.5, the 
limit of the Steidel et al. image, there are no luminous candidates within ~ 30" of the QSO. 

Based upon the residuals following the point-spread function subtraction of the QSO in both the 
WFPC and the ground-based images, there is no evidence for luminous objects directly in front of the QSO 
("zero-impact" absorbers). However, dwarf galaxies of roughly < 0.01 L* at zero impact cannot be ruled 



'^Depending upon the Population III IMF, some pockets of the IGM may have experienced slow metallicity build up due to 
late type stars. However, we find this scenario to be no difTerent in name than if the process took place in "galaxies". 



- 20 - 



out. In general, it seems unlikely that two absorbing galaxies, separated by such a large redshift interval, 
would be aligned with the QSO on the sky. According to the work of Bowen et al. (1997), dwarf spheroid 
galaxies similar to Leo I are not massive enough to have halos that can contribute significantly to the metal 
line absorption cross section of QSO absorbers seen at high redshift. But the situation is not so clear overall, 
given the recent discovery by BBLD of the saturated Mgii doublet associated with the z = 0.072 dwarf 
galaxy at impact parameter 2.1hj^ kpc. The emission line properties of this dwarf ( ^teidel et al. 1993| ) 
suggest that star formation in these objects may directly govern their gas cross section. If so, perhaps active 



star forming dwarf galaxies could contribute to the overall metal line absorption cross section (cf. York 
et al. 1986| ; |Yanny fc York 1992| ). Naively, one would then expect that the abundance pattern arising from a 
bursting dwarf would be a-group enhanced. Also, it is likely that UV ionizing flux from the newly formed 
O and B stars would contribute to the ionization conditions in the absorbers, which is not what we find. 

If we are to assume that these two absorbers are associated with galaxies of some type, and if we 
accept the lack of evidence for HSB candidates in the PKS 0454+039 field, then we must explore the idea 
that low surface brightness (LSB) galaxies (cf. Bothun et al. 1997) could be giving rise to the absorbing 



gas. Particularly, we are lead to consider the class of galaxies called "giant LSB galaxies" (Bprayberry et al 
1993| ; [Sprayberry et al. 1995| ). 



Great progress in our knowledge of LSB galaxies, their number density, sizes, metallicities, and 



luminosity function in the local universe has been made over the past few years (de Blok 1997; Quillcn 



fe Pickering 1997a; Dalcanton et al. 1997; Sprayberry et al. 1997; de Jong 1996| ; ^prayberry et al" 
1995; McGaugh 1994). Dalcanton et al. (1997) find that LSB galaxies have a space density of at least 
n = 0.03 galaxies h^^ Mpc~'^ and outnumber comparable HSB galaxies by factors of 2 or more. LSB 
galaxies are a non-negligible component of the local universe baryonic mass. Sprayberry et al. (1995) find 
that LSB giants have larger disk scale lengths than HSB galaxies of comparable total luminosity. In the 
cases of F568-6 and UGC 6614, the luminous spiral arms extend to 80h^^ and 50h^^ kpc, respectively 



( Quillcn fc Pickering 1997a ). The extent of the Hi disks for the general population of LSB galaxies is seen 
to be roughly 2.5 times that of their D25, the diameters of their /is = 25 mag arcsec"^ isophotes (van dei 
Hulst et al. 1993). If this scaling holds for F568-6, then it may have an Hi disk of ~ 200/1"^ kpc. In the 



case of the LSB galaxy 1226 + 0105, the Hi disk may extend more than four times its D25 (Sprayberry et al 
1995| ). 



Whatever structures give rise to these two absorbers, the Hi gas must have a velocity dispersion of no 
more than ~ 30 km s^^, as dictated by the inferred upper limit on the b parameter of the z — 0.6428 cloud. 
The constraint is even as low as ^ 15 km s^-'^ based upon the z = 0.9315 cloud. The best values of the 
cloud b parameters are ~ 12 km s^^ and ^ 9 km s~^, respectively. These b parameters fall below the lower 
cut offs in the overall Lya forest b distribution ( Kim et al. 1997 ; Lu et al. 1996| ; Hu et al. 1995 ). Again, 
this suggest that these absorber possibly arise in a minority sub-class of the overall Lya cloud population. 
In their study of the giant Hi disks of F568-6 and UGC 6614, Quillen and Pickering (1997a) found that 
the Hi showed small velocity dispersions of 10-30 km s~^ and 10-20 km s""'^, respectively, as compared to 
60-90 km s^^ measured for local HSB spirals (Vogel et al. 198? ; Canzian et al. 1993). These dispersions 



were measured among the spiral arms; it could be that the extended outer disks are even more quiescent. 
Even so, these values are consistent with the allowed ranges of the two absorbers. If a sub-population of the 
Lya forest is arising in LSB galaxies, they may be characterized by having b parameters scattered about 
the low end of the distribution. 

The metallicities of several LSB galaxies have been measured using Hii regions. LSB galaxies with 
relatively smaller disk scale lengths are found to have [Z/Zq] < —0.3, and are therefore metal poor 
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(McGaugh 1994). However, McGaugh found near-solar and super-solar metallicities for UGC 5709 and 
F568-6, respectively. These two galaxies have large disk scale lengths, and classify as giant LSB galaxies.]^ 
Pickering and Impey (1995; C. Inipey, private communication) have also found other giant LSB galaxies 
have metallicities that scatter around solar. These giant galaxies are found to have stellar surface densities 
at least on the same order as their gas densities, which leads Pickering and Impey to suggest that these 
galaxies have been forming stars slowly. 

We find these facts to be quite interesting in light of the two high metallicity Mgii systems we have 
found. It is well established that Mgii absorption with > 0.3 A selects the population of HSB galaxies 



(Churchill, Steidel, fc Vogt 1996; 3teidel 1995; Steidel, Dickinson, fc Persson 1994; Bergeron fc Boisse 



1991). These galaxies appear to be "normal" in their morphologies and to have luminosities greater than 



0.05 L^. By comparison, the total luminosities of giant LSB galaxies scatter about L* (cf. Sprayberry 



et al. 1995). Though LSB galaxies have low luminosity densities, their disks are proportionally larger. 



giving them total luminosities on par with HSB galaxies. Thus, it seems reasonable that LSB galaxies could 
be part of a more general Mgii absorption selected galaxy population, where the LSB galaxies are selected 
by the smaller Mgii equivalent widths. 

In a recent survey to a limiting rest-frame equivalent width of 0.02 A, Churchill et al. (1997) found 
that Mgii absorbers with Wr(A2796) < 0.3 A (hereafter called "weak Mgii absorbers") account for ~ 65% 



of all Mgii absorbers (also see Churchill 1997b ) Nothing is yet known about the type of luminous object they 
select; none have luminous candidates to < 0.06 L*j^ [assuming the Freeman (1970) surface brightness] in the 
survey of Mgii absorbers by Steidel et al. ( p. Steidel, private communication ). These systems also exhibit 
Fell absorption; for the sample, (logiV(Feii)/A^(Mgii)) = -0.3 ± 0.4 (i.e. they may have [Fe/H] > -1). 
Without information on their Hi absorption, we can only speculate that some fraction of the weak Mgii 
systems have ionization and chemical conditions similar to the two absorbers studied in this work. 

If weak Mgii absorbers are selecting out a "missing" part of the Mgii absorption selected galaxy 
population, LSB galaxies are a logical candidate for this missing portion, particularly the class of giant 
LSB galaxies, or "Malin-cousins" , as designated by Sprayberry et al. (1993, 1995). These galaxies are 
disk galaxies, and these disks are observed to have a lower neutral Hi surface density than HSB galaxies 
(Bothun et al. 1997). As a result, the giant LSB galaxy disks have relatively quiescent stellar evolution. In 
fact, the general population of LSB galaxies show a trend of increasing red color with increasing disk scale 



length (Sprayberry et al. 1995). Quillen and Pickering (1997b) reported extremely red colors {R — H = 2.2 
and B — H = 3.5 — 4.2) for the two giant LSB galaxies UGC 6614 and F568-6, which suggest that they 
have a dominating old component in their stellar populations. These galaxies provide the precise type of 
environment in which there has been ample time for iron-group enhancement and metallicity build up in 
the gas phase of their disks. Further, the quiescent nature of their disks leads us to suggest that we should 
not expect to see the complex velocity structures seen in the majority of the stronger Mgll absorption 
profiles ( phurchill 1997a ; Churchill et al. 1998). Indeed, the low Hi surface density of these galaxy disks 
should result in weaker Mgii absorption because there is less of the neutral hydrogen shielding required for 
Mgii to survive. Also, the quiescent nature of the interaction between the gas and stars in these galaxies 
suggests that the gas is not being stirred up, which would generate erratic gas kinematics. It may be that 
such processes facilitate the generation of a high ionization layer around the disk, as seen in the Galaxy 



*The giant, or large scale length, LSB galaxies are defined by Sprayberry et al. (1995) to have fig (0) + 5 log > 27, where 
^J.BiO) is the central surface brightness in the B band and is the scale length in h~^hpc. The high metallicity LSB galaxies 
are seen to have > kpc, which constitute ~ 1/3 of the known giant LSB galaxies defined in Sprayberry et al. 
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(cf. Ravage fc Scmbach 1996| ). The small b parameters inferred for the Hi and the apparent lack of Civ 
absorption in the two weak absorbers are consistent with a quiescent disk. 

If weak Mgii absorbers are selecting giant LSB galaxies, then these absorbers provide a potential 
probe of the number density of these massive galaxies. LSB galaxy disks may grow from isolated l-2a 
peaks in the initial density fluctuation spectrum and may trace low density extended dark matter halos in 
a relatively unbiased way (Bothun et al. 1997). They also would provide a powerful probe of the chemical 
enrichment history of LSB galaxies, which appear to evolve at a significantly slower rate and may produce 
stars via conventional pathways [such as not within molecular clouds (cf. Bothun et al. 1997)]. 

To date, there is not enough known about the number density and gaseous cross sections of the class 
of giant LSB galaxies to compare their dN/dz directly to that of the weak Mgii systems, or to place 
meaningful limits on their number evolution if we assume they are selected by weak Mgii absorption 



(Bothun et al. 1997; C. Impey, private communication). Roughly, the overall population of LSB galaxies 



appear to follow a trend such that those with larger disk scale lengths are observed to have smaller central 
surface brightness ( Sprayberry et al. 1995| ). Following this relation, there is a significant gap between 



Malin 1 and the remaining sub-population of giant LSB galaxies. Malin 1 has a disk scale length a factor 
of 30 times greater than that of F568-6 and a central surface brightness a factor of 100 less than F568-6. 
Does the gap in this parameter space reflect a true break, suggesting that Malin 1 is a rare galaxy type? 
Or, is the gap an artifact of selection effects? As Sprayberry et al. point out, it is important to explore this 
parameter space in order to determine the size and number density distributions of giant LSB galaxies. If 
LSB galaxies are considered to be a natural extension of the HSB galaxy luminosity function, and their disk 
scale lengths and central surface brightnesses exhibit similar behaviors to those found for HSB galaxies, 
then the region of central surface brightness - scale length parameter space giving rise to Malin-type LSB 



galaxies is continuously populated (3. Linder, private communication) 



We can only tentatively suppose that weak Mgii systems with accompanying Fell absorption may be 
selecting the class of giant LSB galaxies ( Sprayberry et al. 1993 ) out to z ~ 1. Impey & Bothun (1989), 



upon reexamining the selection effects and assumptions that go into the calculations of galaxy cross sections 
from QSO absorbers, found that LSB galaxies are expected to dominate the absorption cross section. If 
we assume a non-evolving absorber cross section, noUo, we can write the projected radial extent of QSO 
absorbers as i? ~ 7.5{N/noy^^h^^ kpc, where we have integrated over the redshift interval 0.42 < z < 1.18, 
and where N is the number of observed absorbers in the redshift interval. If we are to claim that the two 
Mgii systems arise from the general population of LSB galaxies, we obtain i? ~ 60 kpc, where we have used 
the space density, n = 0.03 galaxies /lyg Mpc^'^, found by Dalcanton et al. (1997). Supposing that giant 
LSB galaxies comprise 1% of the LSB population, the size of the absorbers is inferred to be i? ~ 200 kpc, 
which compares to the Hi sizes of F568-6, UGC 6614, and UGC 5709 (the "Mahn cousins"). 

The success rate of two weak Mgii absorbers out of 28 Lya clouds leads us to tentatively suggest that 
5-10% of the so-called Lya clouds in the forest at (z) ^ 0.7 will have detectable Mgii absorption to our 
level of sensitivity. Using the dN/dz for the Lya forest from the Quasar Absorption Line Key Project 



( lannuzi et al. 1998 ), this corresponds to non-evolving population with dN/dz ~ 1.5 — 3. Interestingly, 
this is not inconsistent with the number density of dN/dz ~ 1.74 for Mgii absorbers with Wi- < 0.3 A, 
as reported by Churchill et al. (1997). As with the Lya absorbers at z < 1, the weak Mgii systems are 
consistent with a non-evolving population. It could be that our search through the PKS 0454+039 forest 
has picked up the same population of absorbers selected by a fair fraction of the weak Mgii survey, whatever 
that population may be. Fell absorption is present in many of these systems, suggesting that many of 
these absorbers may have chemical and photoionization conditions similar to the two absorbers along the 
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PKS 0454+039 sight line (i.e. they may arise in similar environments with similar evolutionary histories). 



In contrast to post-starbursting dwarfs, which have not formally been ruled out if they both happen 
to be tightly aligned on the sky with the QSO, the few giant LSB galaxies known to date have colors 
consistent with a population of late-type stars ( puillen fc Pickering 1997b| ). This is consistent with the 
inferred ionization conditions of the two absorbers, since early-type stars are ruled out as sources of UV 



flux. Most giant LSB galaxies have active galactic nuclei [narrow emission lines ( Sprayberry et al. 1995 



Bothun et al. 1997)], and this could very well be the case here. There is a point-like/stellar object (#7 in 
LBBD) with impact parameter ~ 45 or 60h^^ (assuming z = 0.6 or z = 0.9.) that has a one-sided spiral 
arm like structure. Though many damped Lya absorbers are seen to be LSB galaxies and low luminosity 
dwarfs ( [Cohen et al. 1996| ; |Le Brun et al. 1997| ; [Steidel et al. 1997] ; [Meyer fc York 1992| ), this is not a 
necessary condition for the absorbers to be LSB galaxies; the high redshift LSB absorbing galaxies studied 
so far have been selected by their Hi absorption and not by their Mgii absorption. Weak Mgii absorption 
that is accompanied by Fell absorption of comparable strength may be selecting a well-defined population 
of luminous objects. These objects probably do not include the damped Lya systems, but probably do 
include sub-Lyman limit systems with high metallicity. This is in contrast to the findings of BBLD, who 
find that strong Mgil absorption accompanied by strong Fell absorption selects damped Lya systems of 
low metallicity. 

The weak Mgii absorbers may in fact be selecting LSB galaxies, given that the weak Mgii absorbers 



found by ChurchiU et al. (1997) do not have HSB candidates to < 0.06 L*j^ ( [Steidel 1995[ ; p. Stcidcl, private 
communication). Very deep imaging and faint object spectroscopy will be required if we are to identify 
the luminous objects selected by weak Mgii absorption. These objects may represent a significant fraction 
of the galaxy population of the universe (and therefore dark matter content). Thus, understanding their 
statistical properties is important for theories of structure formation and galactic evolution. 

In the future, a detailed comparison of the relative abundances of a-group elements (O, Ne, Mg, Si, 
S, Ca) to Fe-peak elements (Cr, Ni, Fe, Zn) in "Lya forest clouds" holds the promise of revealing their 
various origins and formation epochs. The rate of Type la SNe appears to be very low for the first Gyr in 
the history of a Milky- Way like galaxy ( Truran fc Burkcrt 1993 ; Smecker fc Wyse 1992| ). If the delay is 
similar for the onset of Type la SNe in other galaxies as well, then it could be that [Mg/Fe] abundance ratio 
tests are confined to low redshifts. As such, it is expected that a-group enhanced abundance ratios should 
be almost exclusively seen at z > 1.5 (Timmes et al. 1995). For go — 0.5 and A = 0, if a galaxy formed at 
z > 4, then Type la SNe may not begin to contribute to its chemical enrichment until z = 1.5 and may not 
influence the abundance patterns to a detectable level until z '-^ 1 (1 Gyr later). For qo ~ 0.1 and A = 0, 
the Fe-group enrichment would likely be seen no earlier than z ~ 1.5. 



8. Summary 

We searched for Mgii AA2976, 2803 absorption in 28 Lya forest absorbers along the PKS 0454+039 



line of sight. The spectrum studied for the metal line absorption was an optical HIRES spectrum (Churchill 



1997a[ ; ChurchiU et al. 1998). The Lya line fist were taken from the UV G190H and the G270H FOS/HST 



spectra of BBLD. The redshift range was 0.4163 < z(A2796) < 1.1871. The doublets were identified and 
confirmed using the techniques described in Schneider et al. (1993) and in Churchill ct al. (1997). We found 
two Mgii absorbers, one at z = 0.6428, and one at z = 0.9315. Both these systems exhibit Fell absorption. 



We carefully searched the FOS/HST spectrum for the expected metal-line transitions (cf. Hellsten et al. 
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1997) that were covered in the forest. There are currently no other high spectral resolution studies of metals 
in the z < 1 Lya forest, and it also appears that there are no counter parts to these two systems at higher 
redshifts. 

In Figure |l], we present the 5a rest-frame observed equivalent width detection limit of the 
Mgii A2796 transition as a function of redshift. The detection limit rang ed from 0.007 to 0.020 A, 
except for 2:(A2796) < 0.4662, where it ranges from 0.020 to 0.035 A. The 5(7 mean upper limit is 
logiV(Mgii) ^ 11.3 cm~^ for clouds with 0.1 < Wr(Lya) < 1.6 A, which corresponds to neutral hydrogen 
column densities over the range 13.5 < log7V(Hi) < 18.5 cm~^. 

We studied the two discovered absorbers in some detail. Voigt profile fitting was performed to 
obtain the column densities and Doppler b parameters of the Mgii and Fell. Since the Mgii and Fell are 
unresolved, we performed Monte Carlo modeling of the Voigt profile fitting in order to best constrain the 



uncertainties in the Voigt profile parameters (Churchill 1997a; this work). We then used CLOUDY (Ferland 



1996| ) to model the ionization and chemical conditions of the two absorbers. CLOUDY was used in its 



optimize mode, in which the residuals between the model and the measured Mgii and Fell column densities 
were minimized. The fixed quantities for each cloud, which constitute the grid parameters, were the UV 
flux intensity and continuum shape, the metal abundance pattern, and the "observed" neutral hydrogen 
column density. The optimized output quantities were the total hydrogen density, nn, and a scaling factor 
for the metallicity, Zgcaic- We constrained the model clouds using the Hi column density and the parameter 
/ = 6turb/^tot, where bt ui-b is the non thermal contribution to ^tot- 

The UV continua considered were a Haardt & Madau (1996) UV background (UVB) spectrum, several 



Kurucz (1991) Atlas stellar models, late-type galaxy models ( Worthey 1994 ), and a star forming galaxy 



model ( iBruzual fc Chariot 1993 ). We investigated three abundance patterns: solar, Hii depletion, and ISM 



depletion (the latter being a-group enhanced). The Hii and ISM CLOUDY models included grain physics. 
The only UV ionizing scenario that yielded model clouds that were consistent both with astrophysical 
constraints (numbers of stars, etc.) and with constraints imposed by the data was the Haardt & Madau 
UVB. We conclude that the absorbers are photoionized by the UVB and not by stellar radiation. Neither 
absorber is consistent with having an a-group enhanced abundance pattern. 



As described in |6.l| , the z — 0.6428 absorber may have a near-solar or super-solar [Fe/H]. For the solar 
abundance pattern, the model cloud has 16.3 < logiV(Hi) < 16.8 cm^^, a 6turb/^tot of 0.85 < / < 0.93, 
metalhcity —0.2 > [Z/Zq] > —0.7, and density 0.01 < nn < 0.02 cm~^. For the Hii abundance pattern, 
the absorber has 16.7 < logA''(Hl) < 17.2 cm"^, fcturb/^tot in the range 0.90 < / < 0.95, metallicity 
+0.4 > [Z/Zjjjj] > 0.0, and density rtn — 0.008 cm^'^. If this cloud is relatively free of dust depletion, so 
that the abundance pattern is close to solar, then the cloud has [Fe/H] > —1. If the gas-phase abundance 
follows that of depleted clouds in our Galaxy, then the cloud could have [Fe/H] > 0. 



As described in §6.2, the z — 0.9315 model cloud appears to have a super-solar gas-phase [Fe/H]. 
For the solar abundance pattern, the absorber has 15.8 < logiV(Hi) < 16.3 cm^^, a 5turb/&tot of 
0.62 < / < 0.94, metallicity +0.7 > [Z/Zq] > +0.1, and density 0.2 < nn < 0.4 cm^^. The inferred [Fe/H] 
is even greater for the Hii abundance pattern model. For the Hii abundance pattern, the absorber has 
16.0 < log7V(Hi) < 16.5 cm-2, 6turb/&tot in the range 0.75 < / < 0.97, metallicity +1.6 > [^/%ii] > +0.7, 
and density nn ~ 0.1 cm"'^. This is a metallicity enhancement of five to 40 times over the typical values 



seen in Galactic Hii regions (Baldwin et al. 1991; Rubin et al. 1991; Osterbrock et al. 1992). No matter the 



abundance pattern, this cloud has [Fe/H] > 0. 

In §1^, we discussed the possibility that these two absorbers could arise in giant LSB galaxies. These 
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galaxies are seen to have metallicities that scatter about solar ( Pickering fc Impey 1995| ) and large extended 
disks ( Quillen fc Pickering 1997a ). We tentatively suggest that 5% (at most 10%) of z < 1 "Lya forest 
clouds" with 0.1 < Wr(Lya) < 1-6 A will exhibit Mgii absorption to a 5a detection hmit of 0.02 A. The 
sub-sample of these systems that also exhibit comparable Fell absorption may have iron-group enhanced 
metallicities with [Fe/H] > — 1, and may be selecting giant LSB galaxies at high redshifts. 



This work has been supported in part by the National Science Foundation grant AST-9617185 at 
Penn State. CWC acknowledges support through the Eberly School of Science Distinguished Postdoctoral 
Fellowship at The Pennsylvania State University. Thanks to: P. Boisse for providing the FOS spectrum 
prior to publication; J. Charlton for assistance with the Voigt profile fitting simulations and CLOUDY 
modeling; G. Ferland for making CLOUDY a public tool, U. Hellsten for providing the Haardt fc Madau 
input spectra for CLOUDY; S. Linder for generating plots of hb{0) versus log a of LSB galaxies for our 
inspection and for discussions about the observed properties of LSB galaxies; and C. Steidel for sharing 
unpublished data on the PKS 0454+039 field. It is a pleasure to acknowledge J. Bergeron, M. Bershady, 
P. Boisse, J. Charlton, R. Dave, U. Hellsten, C. Impey, J. Lauroesch, P. Petitjean, D. Schneider, and 
M. Shetrone for stimulating discussions and/or comments. Special thanks to S. Vogt for HIRES and for 
providing the opportunity to use it. We thank C. Impey, the referee, for valuable comments that improved 
the quality of this manuscript. For J. L. N. 



APPENDIX 

A. Searching the FOS Spectrum 

It is important to thoroughly check for the presence of absorption from from both low and high 
ionization species in the FOS spectrum. Here, we present our search for absorption in the FOS spectrum 
from both the z ~ 0.6428 and z = 0.9315 systems. We then obtained the limiting column densities (to 
3.5a) using curve of growth analysis, where we have assumed both pure thermal and pure turbulent scaling 
to the measured Mgii b parameters. These limits may be useful for further constraining the chemical and 
ionization conditions of the absorber, even if the sensitivity level in the FOS spectrum is not very high. 
The continuum fit has been adopted from the work of BLLD. We have used interactive software of our own 
that employed the line detection algorithm described in Schneider et al. (1993). Since these transitions are 
embedded in the Lya forest, it is very difhcult to make identifications at a high confidence level. Selected 
results are tabulated in Table ^. 



A.l. The z = 0.6428 Absorber 

A 3.5(7 observed equivalent width detection limit, which corresponds to approximately 0.25 A. Based 
upon preliminary CLOUDY models of this system, we have used the Line Observability Index (LOX) 
of Hellsten et al. (1997) to place limits on those transitions that would most likely be detectable in the 
spectrum. Roughly in decreasing order of their LOX, the covered transitions include Civ AA1548, 1550, 
Sim A1206, Siiv AA1394, 1402, Cii A1335 and A1036, Siv A1063, Siii A1190, Ovi AA1032, 1038, Sin A1260, 
and Sii A1259. 

Both the expected positions of the Sim A1206 and Siv A1063 transitions are coincident with the wings 
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of higher order Lyman scries Unes from the strong z = 1.1538 Mgii absorption system. There is no detected 
Slll A1190, which suggests that the Siv is Ukely not present in absorption. Both the Sill A1260 and Sll 
A1259 transition are not detected, nor are the Sill A1190, 1304, and 1527 transitions. Thus, it appears that 
Sill, Sim, Sii, and Siii are not detected to a 3.5cr observed Umit of 0.25 A, which corresponds to the upper 
hmits log TV ^ 15.4, 15.3, 18.7, and 16.1 cm~^, respectively. Additionally, neither the Alii A1671 nor the 
Aim A1855 transitions are detected, giving limits logiV ~ 14.3 and 14.5 cm~^, respectively. The QSO flux 
is extinct due to the Lyman limit break of the z = 0.8596 damped system at the expected positions of Cii 
A1036 and the Ovi doublet. However, the Cii A1335 transition is not detected, which provides a column 
density limit of log TV ~ 15.2 cm~^ for Cii. We also obtained (not very useful) upper limits on Nil A1084 
and Fciii A1122, giving logiV 16.7 and 17.6 cm-^. 

Determining the limits on the Civ and Siiv doublets is less clear. The CiV A1548 transition, predicted 
to be observed at 2543.37 A, lies precisely at the position of a ^; = 1.0922 Lya line with observed equivalent 
width 0.64 A, as identified by BBLD. There are no other transitions, including Ly/3, that corroborate the 
Lya absorber. We have also measured the wavelength and observed cqiiivalcnt width of the feature and 
obtain Aobs = 2543.27 ± 0.14 A and Wobs = 0.57 ± 0.09 A. There is an unresolved 2.5c7 absorption feature 
0.51 A redward (+90 km s~^) of the expected Civ A1550 transition with Aobs = 2548.11 ± 0.18 A and 
W^obs = 0.22 ± 0.08 A. Though the doublet ratio is consistent with physically allowed values, there are at 
least two reasons to not adopt this feature as the Civ A1550 line a.t z — 0.6428. First, the detection falls 
below a 3.5ct significance when the continuum fit is adjusted downward by only 0.5{S/N)~^, where S/N is 
the signal to noise in the continuum of the normalized spectrum. The continuum placement is not obvious 
to this level of accuracy over the region 2540 — 2558 A, which brackets the putative CiV line. Second, the 
measured Aobs is approximately a 3a \ difference from the predicted Aobs- 

The region of the FOS spectrum where the Siiv doublet transitions are predicted to fall (2289.66 and 
2304.47 A, respectively) exhibits what appears to be a low level fiux depression. The only feature identified 
by BBLD is a broad and asymmetric z = 0.8816 Lya line. The Siiv A1393 transition could be a blend in 
the red wing of this Lya line; the optimized equivalent width spectrum (following Schneider et al. 1993) 
shows that there are two minima, suggestive of a two component blend with Aobs = 2286.4 and 2289.5 A. 
The latter is consistent with the predicted location of the A1393 transition. Additionally, the optimized 
equivalent width spectrum has a single minimum at 2304.7 A, consistent with the predicted location of the 
A1402 transition. These minima have ^ 3.5cr detection significance levels. However, the continuum fit is 
also somewhat uncertain over this region of the spectrum, so that it is difficult to assess if the Siiv doublet 
has actually been observed in absorption. If the SiiV is present, then the CiV is most certainly present. 

All things considered, we adopt the position that the CiV doublet has not been detected in the FOS 
spectrum and assign the 3.5cr observed equivalent width limit of 0.25 A to the Civ A1550 transition, the 
A1548 transition being compromised by the z = 1.0922 Lya line. This gives an upper limit column density 
of log A'' ~ 14.7 cm^^ for CiV. We also adopt the position that the SiiV doublet has not been detected, 
which gives the limiting column density log ~ 15.7 cm~^. However, these are not secure conclusions. If the 
absorber is multiphased, then our analysis is affected only in the sense that our single-phase cloud model 
must be valid only for the cooler phase with the higher density. If so, then the estimated Hi column density 
is an upper limit in this phase of the cloud, since some fraction of it must then arise in the hotter, lower 
density phase. 
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A.2. The z = 0.9315 Absorber 

Again, the 3.5cr observed equivalent width limit of the spectrum is approximately 0.25 A, and we 
have used the LOX of Hellsten et al. (1997) to place limits on those transitions that would most likely be 
detectable in the spectrum. In approximate order of their LOX, the covered transitions include Cm A977, 
Civ AA1548,1550, Sim A1206, Siiv AA1394, 1403, Nm A990, Cii A1335 and A1036, Sm A1012 and A1190, 
Siv A1063, Ovi AA1032, 1038, and Siii A1260, among others. 

None of these transitions were detected to the 3.5cr level. For the clear non detections, the 
corresponding upper limit column densities are logA^(Ciii) ^ 16.3 cm~^, logiV(Siiv) ~ 16.1 cni^^, 
logA^(Niii) - 17.9 cm~^ logiV(Siii) ~ 16.2 cm^^^ logiV(Siii) ~ 15.6 cm"^. For a few of the 
transitions, the limits are affected by the presence of other absorption features in the spectrum. 

The continuum level fit in the region where the SiV A1063 transition is predicted is somewhat uncertain, 
but if the BBLD fit is adopted, then the not very stringent upper limit is logA^(Siv) ^ 17.9 cm~^. The 
Civ doublet, if present, would be blended in both transitions. The A1548 component is coincident with the 
blue wing of the Fell A1608 transition of the z = 0.8598 damped Lya absorber and the A1550 component 



is coincident with the blue wing of the Mgii A2796 transition of the z — 0.0714 dwarf galaxy (see Bteidel 



et al. 1993; Lc Brun ct al. 1997). Thus, we cannot place useful limits on the Civ column density. 



The expected position of the Sim A1206 transition is coincident with a z = 0.8598 Ovi A1038 line from 
the damped Lya absorber. Otherwise the limit would be logA^(Siiii) ~ 15.5 cm^^. The A1206 transition 
is the only transition with which limits on Sim could have been placed. There is clearly no absorption 
feature where the Cii A1036 transition is expected, giving the not very stringent upper limit on Cii of 
logiV ~ 16.8 cm~^. The expected position of the Ovi A1038 transition is coincident with a higher order 
Hi line from the strong Mgii absorption system at z = 1.1537. However, a non-stringent upper limit of 
log ^ 18.6 cm~^ is obtained for Ovi via the A1032 transition, which would arise in a featureless region of 
the spectrum. 

Interestingly, there is a formally significant absorption feature at Aobs = 2168.16 A with 
Wobs = 0.37 ± 0.09 that is coincident with Feiii A1122 at z = 0.9315. Assuming a thermal line 
broadening in a single-phase cloud, this implies the very large column density of log ^ 17.8 cm~^. This 
large of a column density is virtually impossible to reconcile with the Fell measurement. There are two 
alternative possible identifications for the putative Feiii line. First, it could be a Lya line at z = 0.7835, 
though there is no corroborating evidence to support this interpretation. Second, it could be Ly/3 at 
z = 1.1137, which would be supported by the presence of a Lya line at Aobs — 2569.58 A. BBLD detected 
an asymmetric line with an extended blue wing at 2569.85 A, which they identified as a blend of Nv A1242 
at z = 1.0678 (associated with a Civ system) with Sill A1193 at z = 1.1536 (from the strong Mgii system). 
The Nv is the weaker transition of the AA1238, 1242 doublet and the Sin is the stronger transition of the 
AA1190, 1193 doublet. The Nv A1238 and Sin A1190 components of these doublets have been identified by 
BBLD in two adjacent and unambiguous absorption features, so the blend is corroborated by the doublet 
counterparts. If both these features have been correctly identified and the 2569.85 A feature is a blend of 
the doublets, then there appears to be an inconsistency in that unphysical doublet ratios are required. The 
equivalent width of the Nv A1238 - Sin A1190 blend has a lower limit 0.2 A greater than the measured 
equivalent width (under the assumption of optically thick Sin and optically thin Nv). Also, the identified 
blend is asymmetric in such a way as to suggest its centroid may in fact be 2569.6 A (the predicted position 
of the Lya line if the line in question at 2168.16 A were Ly/3). If the blend were actually a Lya line, it 
would place the Nv A1238 and/or the Sin A1190 identifications in question because the observed Lya 
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equivalent width would need to account for at least 50% of the absorption. This particular discussion 
serves to illustrate the level of ambiguity that arises when line identifications in the Lya forest are being 
considered. 

We adopt the position that the line observed at 2168.16 A is not the Feiii A1122 transition at 
z — 0.9315. The strongest argument is that the LOX of the Feiii A1122 transition is significantly smaller 
than that of Cm A977, which has an ionization potential a bit higher than Feiii and is therefore likely to 
be present in an environment giving rise to Feiii. If Feiii were present, the Cm transition would have likely 
been easily seen in absorption. 



B. Computing the Number of Stars 

Assuming no interstellar extinction, stars of all spectral types and luminosity classes produce an 
observed i^fi, = Xf\ = 10"'*'^ ergs cm^^ s^^ at A = 5500 A when their apparent visual magnitudes are 



V = (Allen 1981). Integrating over space out to some finite radius, R, which contains a constant density 



of stars with a given spectral type, one obtains 

log[i./,(0.17)]^^ = log TV, -2 log i? + C,, (Bl) 

where R is in kiloparsecs, and [z//i,(0.17)]^ is the total flux from N^, stars at A = 5500 A, which corresponds 
to 0.17 Ryd. The value of C,, given by log3r//j/(0.17) + 21ogd, — 6, is dependent upon the distance, 

[pc], at which the given stellar type has apparent visual magnitudes V — 0. If the stars are not located 
in a constant density sphere centered on the clouds, but instead are all at some distance, R, the only 



modification to eq. |B1| is that each C, is reduced by a geometric factor, log 3 ^ 0.5. 



The C, are dependent upon the stellar luminosity class, I, III, and V. For 05, BO, AO, and GO stars, 
the C, are 

Co5il,y) = (-5.9,-5.9), 
C7bo(I,III,V) = (-5.8,-6.2,-6.6), 

Cao(I,III,V) = (-5.8,-7.9,-8.5), ^ ' 

Cgo(I,III,V) - (-5.8,-8.6,-10.5). 

Note that all supergiant stars, luminosity class I, have roughly identical C*. 

We are most interested in the flux intensity at 1 — 1.2 Ryd, since this is the energy range corresponding 
to the ionization potentials of the observed transitions. To estimate the number of stars required to elevate 
the integrated stellar /galactic flux such that it equals the UVB vf^, at ^ 1.2 Ryd, we need to introduce a 



continuum shape term into eq. |B1]. We define logfc, to be the ratio of the stellar flux at 0.17 Ryd to that 



at 1.2 Ryd, or log [i^/^ (0.17)]jy — log ['^/iy(l-20)]jY • For a given metallicity. Atlas stellar models (Kurucz 



1991) show that the log/c, are not sensitive to surface gravity (luminosity class), but only to effective 



surface temperature (spectral type). Thus, eq. |B1] can be written 



log A^, - 2 log i? = logfc, + log[i^/,(1.20)]uvB - C*, (B3) 

where log [i^/,y(1.20)]^ has been replaced by log [j^/i,(1.20)]uy3 (under the assumption of no shielding or 
extinction of the UVB), which ranges from 10~^-^ < [i//^(1.20)]uyg < 10~^-^ ergs cm~^ s~-^ in the redshift 
range of interest (see Figure ||) . 

The spectral shapes, especially continuum breaks, play a significant role in the number of stars that 
are required to dominate, or match, the UVB flux at ~ 1 Ryd. For O, B, A, and G stars, respectively. 
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the logfc, are roughly 0, 3, 9, and oo. Thus, for a cloud embedded in a i? = 1 kpc volume of space, the 
number of luminosity class V stars required to influence the ionization conditions of the model clouds are 

^ 1,10'*, 10^^, and oo stars, respectively. For supergaints, the numbers are iV* ^ 1,10'^, 10^, and oo 
stars, respectively. At redshifts of 0.5 < z < 1.0, a single O star could carve out a i? = 1 kpc volume in 
which it could match or exceed the UVB at 1 — 1.2 Ryd. However, it would require 1000 BOI stars, or 10,000 
BOV stars distributed in that same volume. For luminosity class III stars, the numbers are iV* ~ 10^, 10^"'^, 
and oo stars, respectively, where the O type star has been omitted. Because of their sharp Hi break, G, K, 
and M stars would need to have astrophysically unrealistic number densities to contribute to the ionization 
conditions of the absorbers, regardless of their luminosity class. 

It is of interest to know if the implied stellar number densities are consistent with those of the Galaxy. 
In terms of the stellar number density per cubic parsec, n*, eq. |B3] is written, 



logn, +logi? = -9.6 + logfc, +log[i^/,.(1.20)]uvB - C, stars pc'^ (B4) 

where R, as above, is in kpc. In a volume with R = I kpc, the required number density of O stars, whether 
main sequence or supergiant, is hq ^ 10^^ stars pc^"^. For the remaining spectral types, the required 
stellar density is dependent upon their luminosity class. For an AO star, the required number densities 
are ua ~ 1,100, and 100 stars pc~^ for luminosity classes I, HI, and V, respectively. In contrast, the 
observed stellar density of all giants in the solar neighborhood is vanishingly small, and for AOV stars is 



10~3-3 stars pc~3 (Allen 1981) 



Clearly, under the assumption that the UVB is not shielded by layers of neutral hydrogen or 
extinguished due to dust, the type of stellar environment that would be required for the stellar flux to 
dominate over the UVB at 1 Ryd would either need to be uncommonly dense or populated by O stars. 
Lower metallicity stars have smaller breaks and flatter overall continua; lower metallicity stars have smaller 
logfc*. Thus, the solar metallicity stars provide a rough upper limit to the number of required stars. On 
the other hand, interstellar extinction would have the effect of reducing vf^, at 1.2 Ryd for the individual 
stars, which makes the above estimates a lower limit for solar metallicity stars. 
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Table 1. 


Equivalent Width Limits 


in the Lya 


Forest 




Boisso ct al. 1997 




ril 








^abs(Lya) 


W,{Lya) 


log TVgo 


log Wso 


l0giVi5 




log AT 




[A] 


r«w,— 21 
[cm J 


[cm J 


[cm-2] 


[A] 


[cm-2] 


0.4163 


1.03 


14.9 


18.1 


18.5 


0.025 


11.8 


0.4662 


0.37 


14.0 


14.4 


16.7 


0.017 


11.6 


0.5577 


0.36 


14.0 


14.4 


16.5 


0.013 


11.5 


0.5898 


0.63 


14.3 


16.3 


17.8 


0.013 


11.5 


0.6251 


0.35 


13.9 


14.3 


16.4 


0.012 


11.4 


0.6428''^ 


0.33 


13.9 


14.3 


16.4 


0.012 


11.4 


0.6447 


0.75 


14.5 


17.2 


18.1 


0.015 


11.6 


0.6786 


0.28 


13.8 


14.0 


15.4 


0.009 


11.3 


0.6968 


0.19 


13.6 


13.7 


14.2 


0.009 


11.3 


0.7566 


0.29 


13.8 


14.1 


15.5 


0.012 


11.4 


0.7934 


0.11 


13.3 


13.4 


13.5 


0.009 


11.3 


0.7986 


0.11 


13.3 


13.4 


13.5 


0.009 


11.3 


0.8069 


0.39 


14.0 


14.5 


16.9 


0.011 


11.4 


0.8125 


0.35 


13.9 


14.3 


16.4 


0.011 


11.4 


0.8816 


0.53 


14.2 


15.4 


17.6 


0.008 


11.3 


0.9084 


0.59 


14.3 


15.9 


17.8 


0.008 


11.3 


0.9183 


0.78 




1 7 


18.2 


0.011 


11.4 


0.9315^^ 


0.15 


13.5 


13.6 


13.9 


0.008 


11.3 


0.9384 


0.51 


14.2 


15.3 


17.6 


0.008 


11.3 


0.9781 


1.63 


16.5 


18.8 


19.2 


0.008 


11.3 


0.9948 


1.19 


15.3 


18.4 


18.7 


0.008 


11.3 


1.0558 


0.50 


14.1 


15.2 


17.5 


0.008 


11.2 


1.0679 


0.41 


14.0 


14.6 


17.0 


0.008 


11.2 


1.0922 


0.31 


13.9 


14.2 


15.8 


0.007 


11.2 


1.0995 


0.40 


14.0 


14.6 


17.0 


0.007 


11.2 


1.1674 


0.86 


14.6 


17.7 


18.3 


0.007 


11.2 


1.1749 


0.31 


13.9 


14.2 


15.8 


0.007 


11.2 


1.1871 


0.16 


13.5 


13.6 


13.9 


0.007 


11.2 



'^These two systems were detected in Mgii, but are presented here for purpose of 
illustration. 

Note. — Columns 1 and 2 give the redshift and rest-frame Lya equivalent width 
as measured by Boisse et al. (1997), except where noted. Columns 3, 4, and 5 give 
the estimated neutral hydrogen column density for b(Hi) = 80,30, and 15 km s~^, 
respectively. Columns 6 and 7 give the 5a rest-frame Mgii A2796 equivalent width 
limit and Mgii column density assuming linear curve of growth analysis for the 
Wr,ii„i(Mgii). 
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Table 2. Measured Properties of Absorbing Systems 



Transition 


^obs 

[A] 


W, logiV 
[A] [cni-2] 


h 

[km s~^] 






z 


= 0.6428 Absorber 




Hi 1216 
Fell 2383 
Fell 2600 
Mgii 2796 
Mgii 2803 


1997.10 
3914.504 
4271.654 
4593.924 
4605.719 


0.33 ±0.03 
0.031 ±0.012 12.46 ±0.06 
0.029 ±0.004 

0.120 ±0.008 12.74 ±0.02 
0.086 ± 0.004 


4.29 ± 1.25 
5.73 ±0.28 




z 


= 0.9315 Absorber 




Hi 1216 
Fell 2383 
Mgii 2796 
Mgii 2803*^ 


2348.07 
4602.312 
5401.148 
5415.035 


0.15 ±0.04 
0.020 ±0.003 12.24 ±0.09 
0.041 ± 0.002 12.29 ± 0.08 
0.020 ±0.001 


2.28 ± 1.57 
2.19 ±0.46 



'^Bascd upon the Monte Carlo simulations, the equivalent width of 
this transition may actually be 0.023 A. 

Note. — Columns 4 and 5 are based upon Voigt Profile fits to 
the HIRES data. The tabulated column densities and h parameters 
apply to all transitions of a given ion species; however, they are given 
opposite only to the strongest transition for the species. 
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Table 3. FOS Limiting Column Densities of Mgii Absorbers 







Thermal Limits 


Turbulent Limits 


Transition 


Acxp 


b 


logA^ 


b 


logA^ 




[A] 


[km s^-'-] 




[km s-i] 


[cm^^] 




z = 0.6428 Absorber 


Sm A1190 


1955.65 


4.9 


< 16.1 


5.7 


< 15.9 


Sim A1206 


1982.07 


5.3 


< 15.3 


5.7 


< 15.2 


Sin A1260 


2070.66 


5.3 


< 15.4 


5.7 


< 15.3 


Cii A1335 


2194.34 


8.1 


< 15.2 


5.7 


< 16.4 


Siiv AA1394, 1402^ 


2289.70, 2304.51 


5.3 


< 15.6 


5.7 


< 15.2 


Civ AA1548,1550^ 


2543.42, 2574,65 


8.1 


< 14.7 


5.7 


< 15.7 




z = 0.9315 Absorb. 


er 






Cm A977 


1887.11 


3.1 


< 16.3 


2.2 


< 16.4 


Ovi AA1032,1038 


1993.17, 2004.15 


2.7 


< 18.6 


2.2 


< 18.7 


Cii A1036 


2001.68 


3.1 


< 16.8 


2.2 


< 17.0 


Sm A1190, 


2299.29 


2.0 


< 16.2 


2.2 


< 16.2 


Sim A1206 


2330.35 


2.0 


< 15.5 


2.2 


< 15.5 


Sin A1260 


2434.51 


2.0 


< 15.6 


2.2 


< 15.6 


Siiv AA1394, 1403 


2692.04, 2709.45 


2.0 


< 16.1 


2.2 


< 16.1 


Civ AA1548,1550 


2990.34, 2995.31 


2.0 


< 16.7 


2.2 


< 16.8 



These limits on these transitions are not highly certain (see Appendix 
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Table 4. Haardt & Madau UVB: z = 0.6428 Optimized CLOUDY Models 



^scalc 


A(Mgii) 


A (Fell) 


log TT-H 


7V(Hi) 


A^(Hii) 


Temp 


&turb/&tot 


Note'' 




[dex] 


[dex] 


[cni~^] 


[cni~^] 


[cm^^] 


[°K] 










Solar Abundance Pattern; No Grains 


2.0 


0.03 


0.00 


-2.4 


14.50 


15.2 


170 


0.99 


- 1.00 


N 


0.9 


0.04 


0.00 


-1.9 


15.50 


16.3 


720 


0.98 


- 1.00 


N 


0.1 


0.03 


0.00 


-2.0 


16.00 


17.4 


6310 


0.93 


- 0.94 


N 


-0.2 


-0.01 


0.00 


-2.0 


16.25 


17.8 


7940 


0.88 


- 0.92 


Y 


-0.4 


0.03 


0.00 


-2.0 


16.50 


18.0 


10000 


0.87 


-0.90 


Y 


-0.6 


0.00 


0.00 


-1.7 


16.75 


18.1 


10200 


0.86 


-0.89 


Y 


-0.8 


0.03 


0.00 


-1.6 


17.00 


18.3 


11200 


0.85 


-0.88 


Y: 


-1.0 


0.00 


0.00 


-1.5 


17.25 


18.4 


11500 


0.84 


-0.88 


N 


Hii Abundance Pattern; Grains 


3.0 


0.00 


0.00 


-2.9 


14.50 


15.3 


90 


0.99 


- 1.00 


N 


2.0 


0.00 


0.01 


-2.4 


15.50 


16.4 


200 


0.99 


- 1.00 


N 


1.5 


0.00 


-0.01 


-2.2 


16.00 


16.9 


490 


0.98 


- 1.00 


N 


1.2 


0.00 


0.00 


-2.1 


16.25 


17.3 


1170 


0.97 


- 1.00 


N 


0.7 


0.01 


0.00 


-2.1 


16.50 


18.0 


4170 


0.95 


- 0.96 


N: 


0.4 


0.00 


0.00 


-2.1 


16.75 


18.3 


6310 


0.92 


-0.94 


Y 


0.2 


0.00 


0.00 


-2.1 


17.00 


18.6 


8130 


0.90 


-0.92 


Y 


-0.1 


0.00 


0.00 


-2.1 


17.25 


18.8 


10000 


0.88 


-0.90 


N: 



^Optimal logarithmic scaling factor applied to all elements heavier than helium for the given 
abundance patterns. 

''Notes: (Y=yes, N=no) indicate whether CLOUDY conditions are consistent with inferred 
/ = ^'turb/^tot and allowed values of iV(Hi) and 5tot shown in Figure ^. A ":" indicates a 
borderline case. 
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Table 5. Haardt & Madau UVB: z = 0.9315 Optimized CLOUDY Models 



^scalc 


A(Mgii) 


A (Fell) 


log riH 


A^(Hi) 


A^(Hii) 


Temp 


^turb /^tot 


Note'' 




[dcx] 


[dex] 


[cm^^] 


[cm^^] 


[cm^^] 


[°K] 










Solar Abundance Pattern; No Grains 


2.0 


0.03 


-0.03 


-1.1 


14.50 


14.5 


100 


0.98 


- 1.00 


N 


0.9 


0.04 


-0.04 


-0.8 


15.50 


15.6 


420 


0.95 


-0.98 


N 


0.7 


0.04 


-0.03 


-0.6 


15.75 


15.9 


810 


0.90 


- 0.96 


Y 


0.3 


0.03 


-0.03 


-0.6 


16.00 


16.3 


3310 


0.50 


- 0.82 


Y 


0.1 


0.03 


-0.04 


-0.4 


16.25 


16.6 


4790 


0.00 


- 0.73 


Y 


-0.1 


0.03 


-0.03 


-0.2 


16.50 


16.7 


5890 


0.00 


-0.63 


N 


-0.7 


0.03 


-0.03 


-0.2 


17.00 


17.3 


8910 


0.00 


-0.38 


N 


Hii Abundance Pattern; Grains 


2.9 


0.00 


0.00 


-1.8 


14.50 


14.8 


80 


0.98 


- 1.00 


N 


2.0 


0.00 


0.00 


-1.2 


15.50 


15.6 


110 


0.98 


- 1.00 


N 


1.8 


0.00 


0.00 


-1.1 


15.75 


15.8 


130 


0.97 


- 1.00 


N 


1.6 


0.00 


0.00 


-0.9 


16.00 


16.0 


200 


0.96 


- 0.99 


Y 


1.3 


0.01 


0.01 


-0.8 


16.25 


16.3 


440 


0.94 


- 0.97 


Y 


0.7 


0.01 


0.00 


-0.8 


16.50 


17.1 


3890 


0.35 


- 0.79 


Y 


0.2 


0.00 


0.00 


-0.6 


17.00 


17.6 


6920 


0.00 


- 0.55 


N 



^Optimal logarithmic scaling factor applied to all elements heavier than helium for the given 
abundance patterns. 

''Notes: (Y=yes, N=no) indicate whether CLOUDY conditions are consistent with inferred 
/ — btnrh/btot and allowed values of Af(Hi) and 6tot shown in Figure ^. 
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.4 .6 .8 1 1.2 

Redshift 



Fig. 1. — The 5ct rest-frame equivalent width detection hmit of the Mgii A2796 transition as a function 
of redshift. The decrease in sensitivity toward lower redshift is due to the HIRES throughput and CCD 
efHciencies. The higher frequency features are due to the blaze efficiency function of the individual echelle 
orders. Gaps in the redshift coverage are appear at 2; 0.83 and above. 
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z = 0.642J 



z = 0.9315 



-H I 1215 1 blend 




-H I 1215 1 




1 









- FOS 






1980 2000 2020 




2320 2340 2360 2380 




3912 3914 3916 

Fe II 2600 




4600 4602 4604 



Fe II 2600 
Lost in the BLOB 
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Mg II 279t 
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■liifvir™vyv+ 
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Wavelength, A 



5412 5414 5416 5418 

Wavelength, A 



Fig. 2. — The FOS/HST (top panels) and HIRES/Keck (metal lines) spectra of the two absorbers. Note 
the very different wavelength scales presented between the FOS and the HIRES spectra. The ticks mark the 
predicted positions of Lya A1215 absorption lines based upon the accurate redshifts of the metal lines. The 
z = 0.6428 Lya line is a member of a blended feature (see Fig. |^). 
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Fig. 3. — The deblending fit for the z = 0.6428 Lya line (second from the left) in the FOS spectrum. The 
blended lines are Ly7 a,t z = 1.1536, Lya at z = 0.6429, and Lya at ^ = 0.6448. The lone feature at 
2004.5 A is Ly6 at ^; = 1.1536, which is the redshift of a strong Mgii absorber. 
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Fig. 4. — The Haardt and Madau UVB normalization and shape for redshifts 1.0 and 0.5, as labeled. The 
ionization potentials for selected ion species are shown as vertical ticks in the range E >1 Ryd. 
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Fig. 5. — The inferred temperature of the z = 0.6428 absorber (upper curves) and of the z = 0.9315 absorber 
(lower curves) as a function of the fraction of the measured Mgli b parameter due to possible turbulence, 
/ = ^turb/^tot- The uncertainties in T, based upon the la uncertainties in the Mgii &tot, are given by the 
dotted curves. The shaded regions give the typical range of T for various objects in the Galactic ISM and 
for Mgii QSO Absorption Line (Mgii/QAL) clouds at z ^ 1. The observed absorbers are consistent with 
Mgii/QAL clouds, Hii Regions, and Warm Hi clouds (see text). 
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Fig. 6. — The z = 0.6428 absorber properties — (a) The range of total Doppler b parameters as a function 
f — foturb/^tot- The dotted Unes are the spread due the la uncertainties in the measured Mgii 6tot- — 
(b) and (c) The range of Hi column densities as a function of /. The thick solid curves are best values of 
A^(Hi) based upon &tot(Mgll) and Wr(Hi). The thin curves give the spread of 7V(Hi) based upon the la- 
uncertainties in 6tot(Mgii) [inner] and Wi(LyQ:) [outer]. The curves that originate in the lower right hand 
corners and rise upward and then to the left are the allowed locus of / for a cloud model with a given iV(Hi). 
Panel (b) is for a solar abundance patterns and (c) is for the Hii abundance pattern. 
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